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Several simple models of strongly correlated bosons on three-dimensional lattices have been shown 
to possess exotic fractionalized Mott insulating phases with a gapless 'photon' excitation. In this 
paper we show how to view the physics of this 'Coulomb' state in terms of the excitations of proxi- 
mate superfluid. We argue for the presence of ordered vortex cores with a broken discrete symmetry 
in the nearby superfluid phase, and that proliferating these degenerate but distinct vortices with 
equal amplitudes produces the Coulomb phase. This provides a simple physical description of the 
origin of the exotic excitations of the Coulomb state. The physical picture is formalized by means of 
a dual description of three-dimensional bosonic systems in terms of fluctuating quantum mechanical 
vortex loops. Such a dual formulation is extensively developed. It is shown how the Coulomb phase 
(as well as various other familiar phases) of three-dimensional bosonic systems may be described in 
this vortex loop theory. For bosons at half-flUing and the closely related system of spin- 1/2 quan- 
tum magnets on a cubic lattice, fractionalized phases as well as bond or 'box' ordered states are 
possible. The latter are analyzed by an extension of techniques previously developed in two spatial 
dimensions. The relation between these 'confining' phases with broken translational symmetry and 
the fractionalized Coulomb phase is exposed. 



I. INTRODUCTION 

Several recent studies have produced a variety of mod- 
els that exhibit quantum number fractionahzation in two 
dimensions and in the absence of magnetic fieldsiSiSi^ 
It has been appreciated for some time that similar frac- 
tionahzation phenomenon can occur also in three dimen- 
sions, and can take novel forms not possible in 2DiSi2i2i^ 
In particular, so called Coulomb phases have recently 
been demonstrated in 3D bosonic models. In these 
Coulomb phases the fractionalized degrees of freedom in- 
teract with a gapless emergent C/(l) gauge field. The 
gauge excitations describe, at low energies and long dis- 
tances, a linear dispersing transverse mode which has 
been dubbed an 'artificial photon'. The adjective 'artifi- 
cial' refers to the fact that this photon really corresponds 
to some collective excitation of the boson system and 
has nothing directly to do with the true electromagnetic 
photon. In particular the 'artificial photon' in the frac- 
tionalized 3D Coulomb phases arises even in microscopic 
models of correlated bosons with purely short-ranged in- 
teractions. The Coulomb phase also contains additional 
gapped excitations that may be identified with magnetic 
monopoles of the internal gauge field. The possibility of 
such dynamical generation of 'light' was noted more than 
two decades ago in Ref. 10. 

In the case of fractionahzation in 2D bosonic systems, 
a very fruitful perspective is obtained by departing from 
a superfluid state and considering nearby Mott insulator 
phases as vortex condensates liitiSiiiii Conventional in- 
sulating states of bosons, with possibly charge or bond 
order, can be understood in these terms. The same per- 
spective also naturally predicts occurrence of a Z2 frac- 



tionalized state when the superfluidity is destroyed by 
the condensation of doubled vortices while single ones 
remain gapped. The unpaired single vortex survives in 
the resulting insulating phase as a gapped excitation that 
retains a Z2 quantum number (only the oddity of vortic- 
ity remains well deflned) — this is the Z2 vortex or vison 
excitation of this fractionalized insulator. 

Turning to 3D systems, we envisage a similar perspec- 
tive on insulating states as vortex condensates, except 
that vortex excitations in the superfluid are now lines in- 
stead of point particles. For example, a Z2 fractionalized 
insulator is obtained by destroying the superfluid order 
by pairing and condensing vortex lines. 

The goal of the present paper is to develop a similar 
intuition for the 3D Coulomb phase (also dubbed U{1) 
phase). In particular, we want to understand the genesis 
of the excitation structure of this phase in an approach 
that departs from the superfluid state. We note that 
the excitation structure of the Coulomb phase consists 
of gapped charge-1/2 chargons that are minimally cou- 
pled to an emergent gapless "photon" mode. In addition, 
there is another gapped excitation that may be viewed 
as a "magnetic monopole" that acts as a source of the 
emergent magnetic flux. The picture that appears from 
the present study is that physics inside vortex cores is 
intimately involved in producing the U{1) phase. We 
thus have a rather unusual but very interesting situation 
when we cannot ignore the microscopies of the vortex 
core when describing the relationship of the superfiuid 
phase to the nearby insulating states. 

The main idea is very simply illustrated by consider- 
ing bosons on a three-dimensional cubic lattice and at an 
average filling of one half boson per site. In some semi- 
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classical description, there are two simple phases that 
might be imagined for the bosons. First, there is a trans- 
lationally invariant superfluid phase. Second, there is a 
Mott insulator in which the bosons preferentially occupy 
sites on one or the other of the two sublattices. Clearly, 
such a Mott insulating ground state is two-fold degen- 
erate corresponding to the two possible states of boson 
charge ordering. The superfluid phase possesses vortices 
which are extended line defects. The superfluid order pa- 
rameter will be suppressed in the core of such a vortex 
line. If the vortex core size is larger than the lattice spac- 
ing (as might well happen) it makes sense to ask about 
what 'phase' is obtained in the core. An immediate guess 
is that the core is simply the checkerboard charge-ordered 
Mott insulating state as illustrated in Fig.^ However the 
checkerboard state breaks a symmetry (that of sublattice 
interchange), and is two-fold degenerate. We are then led 
to conclude that there will be two kinds of vortices de- 
pending on which sublattice has higher boson density in 
the core. 

This simple-minded discussion on the structure of vor- 
tices will be modified once we take quantum fluctuations 
into account more seriously in the vortex core. In partic- 
ular, as the core has finite radius, it essentially behaves 
as a one dimensional system. Quantum fluctuations will 
tend to restore the broken symmetry in the core and pro- 
duce a unique vortex line. However, as the broken sub- 
lattice symmetry is discrete, it is natural to expect that 
the symmetry breaking will be stable to weak quantum 
fluctuations. Thus there will be some range of parame- 
ters in which there will be two distinct kinds of vortex 
lines with symmetry broken cores. In other range of pa- 
rameters a unique vortex with no symmetry breaking in 
the core may well result. The change between the two pa- 
rameter ranges involves no ground state phase transition 
but only a change in the nature of the excited states. 

Now consider disordering the superfluid by proliferat- 
ing and condensing vortex line loops. Clearly, the insu- 
lating phase that obtains will depend on the structure of 
the vortices that are actually proliferating. Consider the 
situation where the vortices have ordered cores. If we 
preferentially condense one or the other kind of vortex, 
we will clearly induce checkerboard charge order in the 
insulator. But what if we condense both species of vor- 
tices with equal amplitude? We argue that the resulting 
phase is the Coulomb fractionalized insulator. We show 
how the excitations of this exotic insulator may be under- 
stood in terms of the various excitations of the superfluid 
when it has ordered vortex cores. 

Consider a vortex line. It is still possible for the core 
order to change from one state to the other somewhere 
along the line as shown in Fig. |21 but there is a gap 
for such a domain wall excitation. We expect this gap 
to remain when such vortices proliferate while retaining 
their order. When both vortex species proliferate, the 
domain wall excitations become point particles which we 
identify with the monopoles of the Coulomb phase. 

The appearance of the artificial photon is more subtle. 



but is also a consequence of proliferating two species of 
vortices. To explain this, it is useful to recall what hap- 
pens in the more familiar situation of particle conden- 
sates. When the particle that is condensing is 'neutral' 
(in the sense of having only short ranged interactions), 
the resulting condensate supports gapless Goldstone ex- 
citations. On the other hand, if it is charged (in the 
sense of having long ranged interactions), there are no 
gapless excitations associated with the condensation (the 
Anderson-Higgs phenomenon). Consider now the dual 
description of bosonic Mott insulators in three spatial di- 
mensions as 'condensates' of extended vortex loops. For 
such a line condensate we may again expect that if the 
vortices have long ranged interactions then the resulting 
condensate supports only gapped excitations. However, 
consider now the situation envisaged above where there 
are two vortex species (distinguished by the ordering in 
the core) which both condense. We may form a 'roton' 
line from these two species by bringing together a vortex 
with one ordered core and an antivortex with the other 
core, and this is pictured in Fig. |21 Such objects have 
no net vorticity, and therefore only short ranged interac- 
tions with one another. When both vortex species con- 
dense such roton loops condense as well. By analogy to 
particle condensates we may expect then to find gapless 
excitations. We will show in this paper that this is indeed 
the case, and that the resulting gapless modes may be 
identified with the photon of the Coulomb phase. More 
precisely, we will argue that the roton lines described 
above may be viewed as the magnetic field lines of the 
emergent deconfined U{1) gauge field that obtains in the 
Coulomb phase. 

It is useful to point a connection with a recent studjii^ 
of so called deconfined critical points in two dimensions 
such as separating a superfluid and a bond density wave 
insulator of bosons at half-filling on a square lattice. In 
2D, a vortex is a point particle, and because of the quan- 
tum tunneling events we expect a unique physical vortex 
in the supcrfiuid phase. In other words, generically the 
two different ordering patterns in the core will tunnel 
into each other at a finite rate in d = 2. However, at the 
critical point at which vortices proliferate, the tunneling 
becomes irrelevant We may then effectively speak of 
two vortex species. The tunneling becomes again rele- 
vant in the vortex condensate, and produces a particular 
insulating state which in this example has bond density 
wave order. However, as elaborated in Ref. |15| there is a 
large intermediate region of length scales in the insulator 
(close to the critical point) which may be described as 
a two-dimensional version of the Coulomb phase. The 
(intermediate) long wavelength physics of this region is 
again described in terms of a gapless excitation which 
may be regarded as the two-dimensional version of the 
photon. In three dimensions, having an ordered vortex 
core is a stable possibility in the superfluid. As men- 
tioned in the previous paragraph, proliferating both vor- 
tex species with equal amplitude then leads to a stable 
insulating phase which has a gapless photon excitation — 
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FIG. 1: View of the cross section of a vortex line in a three- 
dimensional superfluid of bosons at half-filling on a cubic lat- 
tice. Two checkerboard ordering patterns — where the boson 
density is higher on one of the two sublattices — are possible 
inside the vortex core. 




FIG. 2: Vortex lines with different core orders are indicated 
as A (red) and B (light blue). Left panel: A domain wall 
associated with the order in the core. If the superfluid is dis- 
ordered by proliferating both vortex species simultaneously 
while keeping such domain walls gapped, the Coulomb phase 
results. The gapped domain walls survive as the monopoles of 
the Coulomb phase. Right panel: A "roton" loop is formed by 
bringing together a vortex with one ordered core and an an- 
tivortex with the other core. These particular roton loops cor- 
respond to loops of emergent magnetic flux which are strongly 
fluctuating in the Coulomb phase. The artificial photon is as- 
sociated with fluctuations of these loops. 



precisely the Coulomb state at the focus of the present 
study. 

The primary purpose of this paper is to develop the 
physical picture sketched above for the excitations of the 
Coulomb phase in some detail. This goal is strongly 
aided by a formal dual description of three-dimensional 
bosonic systems in terms of extended fluctuating vortex 
loops. This duality transformation — a generalization 
to (3-f 1)D of the familiar boson- vortex duality in two 
spatial dimensions — has been available in the literature 
for many years, and is reviewed in Section llVI Here we 
will use this and further develop it extensively to ob- 
tain a dual description of the phases of various three- 
dimensional bosonic systems which puts meat into the 
physical pictures. In the course of these studies, we de- 
velop a description of bosons at half-filling on a three- 
dimensional cubic lattice that is capable of describing all 
the insulating phases discussed above — the Coulomb 



phase, the checkerboard charge ordered phase, and also 
different bond density wave ordered phases. This gives 
a three-dimensional ^generalization of the corresponding 
development in Ref. 15. Boson systems at half- filling are 
closely related to spin-1/2 quantum magnets. Thus some 
of our results have direct implications for quantum spin 
systems on cubic lattices as well. Indeed we will borrow 
from techniques familiar from studies of quantum anti- 
ferromagnets in two spatial dimensions, and generalize 
them in discussing the structure of the bond ordered in- 
sulating phases of bosons at half-filling 

The paper is organized as follows. In Section m we 
first review a simple bosonic model that explicitly real- 
izes the fractionalized Coulomb phase. This provides a 
good starting point for the discussion of the relationship 
of the U{1) state to other phases. The model has inte- 
ger boson filling and is simple to analyze in the charge 
representation, but the vortex physics is obscured by the 
nontrivial lattice structure. Still, the model motivates an 
effective theory, formulated as a compact C/(l) gauge the- 
ory with two matter fields (chargons) , that properly cap- 
tures the possible phases and is amenable to an analytic 
description in terms of vortices, furnishing the main argu- 
ment for the proposed physical picture. In Sect ion UTTl we 
consider vortex phenomenology in this two chargon the- 
ory. Precise duality analysis is presented in Section llVI 
where we first review the conventional boson-vortex du- 
ality in (3+l)D from a modern perspective, and then 
describe how to include the physics of two vortex species 
in the dual description (the details of the derivation are 
given in Appendix fX)l . In Appendix IdI we consider a 
(2-|-l)D Ising version of the specific boson model in or- 
der to give some intuition in a more simple setting; this 
toy model realizes a topological Z2 phase in two dimen- 
sions via a mechanism that resembles the one proposed 
for the Coulomb phase in three dimensions, with the vor- 
tices replaced by Ising domain walls. In Section Ivl and 
supporting Appendices IbI and ICl we develop a dual de- 
scription for bosons at half-integer filling. An important 
ingredient of this description is monopole Berry phases 
generalizing Haldane's result to three-dimensional spin- 
1/2 quantum magnets on a cubic lattice. We apply this 
in Section IVII where we show how bond density wave 
states can be analyzed as monopole condensates coming 
from the Coulomb phase. 



II. REALIZATION OF THE COULOMB PHASE 
AND PHENOMENOLOGY 

A. Review of boson Hubbard model 

We begin with a review of a specific model that realizes 
Coulomb phase in 3DiS*Si& The model is formulated in 
terms of soft-core bosons (quantum rotors) residing on 
link-centered sites ("islands") of a simple cubic lattice 
as depicted in Fig. 13 The boson islands can also be 
viewed as forming a lattice of corner-sharing octahedra. 
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FIG. 3: Explicit model that realizes Coulomb phase. Boson 
islands are located on the link-centered sites of a simple cubic 
lattice. The islands can be also viewed as forming a network 
of corner-sharing octahedra, and two octahedron units are 
shown. In the Hamiltonian Eq. Q, we stipulate a term UN^ 
which prefers charge neutrality of each octahedron. 

Bosons can hop with amplitude w between neighboring 
islands. Besides the usual on-site repulsion w^, we also 
include repulsive interactions between the bosons that 
favor charge neutrality of each octahedron. The complete 
Hamiltonian is 

H = -w^i^li^v +h.c.) + u^^nf + U^N'^ , 

(W) I r 

where ipj = e**^' creates a boson on a given island and 
ni is the corresponding boson number operator. Each 
octahedron is labeled by the cubic lattice site r at its 
center; the corresponding operator Nr is defined through 

Nr=Y,^l^ (1) 
ler 

where / G r sums over all cubic lattice links emanating 
from r. The total boson number of the system is A^tot = 

We summarize the analysis of the possible phases of 
the above Hamiltonian. When the boson hopping dom- 
inates, w ^ u^,U, the system is a superfluid. When 
the charging energy dominates, u^,U w, the system 
is a conventional Mott insulator. On the other hand, 
as reviewed below, when the charging energies U and 
are varied separately, there is an intermediate regime 
U ^ w ^ ^JuipU such that the system is in a Coulomb 
phased. This phase is a fractionalized insulator with 
the excitation spectrum consisting of gapped charge- 1/2 
chargons, a gapless linearly dispersing photon, and a 
gapped monopole. A schematic phase diagram of our 
model is shown in Fig. ^ 




FIG. 4: Schematic phase diagram of the boson model Eq. Q 
exhibiting the stable phases discussed in the text; note that 
we do not know the details for intermediate coupling strengths 
and whether some other states may intervene. 

The analysis in the limit U ^ w,u^ is similar to that 
in the large U limit of the electronic Hubbard model at 
half-filling. For w — — there is a degenerate mani- 
fold of ground states specified by the requirement Nr = 
for each r, separated by a large charge gap U from the 
nearest sectors. Including the w,u^ terms lifts the de- 
generacy in each such sector, and this is best described 
by deriving the corresponding effective Hamiltonians for 
small perturbing couplings. 

In the ground state sector iV^ = for all r, an elemen- 
tary calculation gives 

( □ 

with ii'ring = 'i'uP' /U . When writing the 'ring exchange' 
term around a given placket, we label a boson island by 
the corresponding link end-points. 

A simple change of variables shows that together 
with the constraint A^^ = can be regarded as the well- 
studied (3-t-l)D compact U{1) gauge theory. Indeed, di- 
vide the underlying cubic lattice of Fig. |31 into A and B 
sublattices. We now define a vector field a,ct = rjrcfir.r+a, 
where rjr = +1 if r g A and r]r = —1 ii r € B, a = x,y^ z; 
we also perform the corresponding transformation to the 
vector field Era = Vr'>T-r,r+a Conjugate to Ura- In the new 
variables 

Nr = rir'^ -E . (3) 

If interpreted as a compact U{1) vector potential 

and Era as the corresponding electric field operator, the 
ground state sector constraint Nr = is simply the Gauss 
law, while the effective Hamiltonian has the declared lat- 
tice gauge theory form 
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In (3+l)D, the compact U{1) gauge theory has two 
distinct phases: For ifping ^ , the gauge theory is con- 
fining, and all excitations carrying non-zero gauge charge 
are confined. Zero gauge charge excitations of course ex- 
ist with a gap oc U; these carry integer physical charge 
in units of elementary boson charge qt- This is the con- 
ventional Mott insulator of our boson model. 

In the opposite regime, ifj-ing ^ u^, the gauge theory 
is in the deconfined Coulomb phase, and has a gapless 
linearly dispersing gauge boson (photon) and a gapped 
topological point defect (monopole) as its distinct excita- 
tions in the Nr — sector. In the charged sectors, objects 
with Nr = I at some site, i.e. physical charge ^ot 
confined and can propagate above a finite gap of order 
U. These charged excitations interact via an emergent 
long-range Coulomb interaction. A detailed description 
of the ground state properties of the Coulomb phase due 
to the gapless photon can be found in Rcf. 8. 

This completes the review part of our discussion, and 
we now focus on the relationship of the U{1) phase to 
the nearby more conventional phases. A two-dimensional 
quantum Ising version of the corner-sharing octahedra 
model is considered in Appendix lUl 



B. Phenomenological gauge theory description 

We want to obtain a description of the Coulomb phase 
that includes the chargon fields; this is required if we 
want to discuss the transition from the Coulomb phase 
to the superfluid as the charge gap collapses to zero. To 
this end, we need to consider charge carrying excitations 
in more detail. In the large U limit, this is done by de- 
riving effective Hamiltonians in the charged sectors, just 
as we did in the ground state sector Eq. Q. For exam- 
ple, to study the motion of a single chargon in a given 
region, we consider the sector with Nr = Srro where tq 
can be anywhere in this region. By inspecting the possi- 
ble moves of the chargon, we find that it can hop only on 
the sites of the same sublattice of the cubic lattice. More 
generally, the total chargon number on a given sublattice 
is conserved for any chargon motion that derives from 
the microscopic boson hopping. This means that there 
are two distinct chargons, and can be also seen directly 
from the Gauss law Eq. An excitation with the gauge 
charge -1-1 carries boson charge +qb/2 if it resides on the 
A sublattice and —qb/'2 if it resides on the B sublattice, 
and in the Coulomb phase the gauge and boson charges 
are both good quantum numbers. 

Instead of working directly with the microscopically 
derived effective Hamiltonian for the chargons, which is 
complicated, we consider the following model gauge the- 
ory 

r rr' 

+ue^E'^ - K^cosiV y. a) . ( 



This has two chargon fields minimally coupled to the 
compact U(l) gauge field. The Hilbert space of the the- 
ory is defined by 

V ■ E ^ nir + n2r ■ (6) 

b\r and b\r both carry gauge charge -1-1, and we take the 
first to carry boson charge +qb/'2. and the second —qb/'^.. 
There is no discrimination between the two chargons, so 
the hopping amplitude and the on-site interaction are 
taken equal for the two species. This theory has the cor- 
rect chargon content and is therefore expected to capture 
the relevant physics of the microscopic model. 

For later convenience, we also write down the corre- 
sponding (3+l)D classical action 

S2ch = [cos(V^(/)i - a^) + cos(Vp02 - a^)] 

-A'^cos(V^at, - V^a^) . (7) 

For simplicity, the action is written in a space-time 
isotropic form, and is characterized by two coupling con- 
stants (3 and K. The former characterizes the relative 
strength of the chargon hopping t vs repulsion C/, while 
the latter characterizes the competition of AT vs u^; . 

As an aside, we note that the above gauge theory 
also appears in a particular "slave boson" treatment of a 
generic bosonic Hamiltonian at integer filling 

i?gonoric = Ub'^nl ~ Wb'^ blbr' H . (8) 

r tt' 

Specifically, consider an analog of CP^ representation for 
the bosons treated as 0(2) rotors 

bl = e*"^" = b\rb2r = e'^"^--"^-) , (9) 

rir — nir — — 7l2r- (10) 

Applying the corresponding slave boson mean field 
scheme, and studying fluctuations about the mean field, 
one obtains the above gauge theory (the details of the 
derivation can be found e.g. in Appendix E of Ref. T^). 
The constraint Eq. (|10|l formally corresponds to the limit 
Mb — !■ cx) in i?2ch, but we can imagine that some coarse- 
grained is finite, and the description in terms of bir 
and b2r as emergent degrees of freedom applies in some 
regime of parameters. 

The phase diagram of the two-chargon gauge the- 
ory H2ch is established by standard argumentsi^ and is 
shown in Fig. |31 The different phases are readily iden- 
tified with the phases of the microscopic boson Hamil- 
tonian Fig. 21 (i) In the regime t U, K ^ ue, the 
gauge theory is in the deconfined Coulomb phase with 
gapped chargons; this is of course the Coulomb phase 
of our microscopic model and has the correct particle 
content by construction, (ii) In the regime t <^ U, 
AT <C Ue, the gauge theory is confining and there are no 
free gauge charges; only gauge-neutral excitations are free 
5)md carry integer boson charge; this phase corresponds to 
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FIG. 5: Schematic phase diagram of the compact U{1) gauge 
theory with two matter fields, Eq. The dashed fine in the 
superfluid phase represents the conjectured vortex core tran- 
sition from a unique vortex to two physicaUy distinct vortices 
(there is no transition in the superfluid ground state across 
this line). 



the Mott insulator, (ill) Finally, in the Higgs phase t ^ U 
both chargon fields condense and the gauge-neutral field 
gj(0i-02) obtains an expectation value; this corresponds 
to the superfluid phase of the bosons. It should be em- 
phasized here that while the model gauge theory repro- 
duces the phase diagram of the physical system, it is only 
in the Coulomb phase that the chargon variables repre- 
sent low-energy degrees of freedom; in the other phases, 
the chargons are confined and are not present as individ- 
ual particles in the spectrum at any energy. 

The phenomenological description allows us in partic- 
ular to discuss the nature of the phase transitions. Thus, 
we expect the conventional insulator to superfluid tran- 
sition to be in the (3-f 1)D XY universality (bosons at 
integer filling) for small K/ue, but it can also become 
first order for larger K/ue- The deconfinement transition 
between the two insulators is similar to that in the pure 
gauge theory and empirically is believed to be first order. 
Finally, in the Coulomb to superfluid transition the gauge 
field can be viewed as noncompact since the monopoles 
remain gapped and do not participate. This transition is 
expected to be fluctuation induced first order lil 

Our treatment so far has been entirely in the boson 
charge degrees of freedom, and to conclude this discus- 
sion we point out some salient properties of the Coulomb 
phase. First, the Coulomb phase features an emergent 
gauge charge conservation law that is exact as long as 
the system remains in this phase. On a formal level, in 
this case, the projective transformation that has led us 
to the U{1) gauge theory is convergent as a series in t/U, 
and the result is a generic C/(l) gauge theory with the 
local constraint that is exact by construction. Second, 
the Coulomb phase and the emergent conservation law 
are stable against adding arbitrary perturbations to the 
microscopic Hamiltonian as long as these perturbations 



are sufficiently small)^iSiLi£ this again follows from the 
local constraint imposed by the large U term, since the 
result is by fiat a gauge theory. This amazing fact holds 
for perturbations that can break all symmetries of the 
microscopic Hamiltonian, and despite the fact that the 
low-energy theory has a gapless photon. This is a hall- 
mark of the topological structure present in the Coulomb 
phase. 

However, the topological character of the Coulomb 
phase and its genesis appear mysterious in the charge 
language. Much further insight is obtained when we con- 
sider vortices in the superfluid phase from the point of 
view of the phenomenological theory. 

A precise topological characterization of the Coulomb 
phase is provided by the existence of the gapped 
monopole excitation. At energy scales well below the 
monopole gap, the physics is that of a non-compact U{\) 
gauge theory. Monopoles remain gapped across the tran- 
sition to the superfluid phase — as we show below this 
leads to unusual vortex physics. 

Specifically, we argue for the presence of two distinct 
stable vortices in the superfluid adjacent to the Coulomb 
phase. We can think of the two vortices as having dis- 
tinguishable order in the core. Physically, the core can 
change from one order to the other somewhere along 
the vortex line, but this domain wall costs an additional 
energy. The stability of two distinct vortices in some 
parameter regime is a consequence of some microscopic 
energetics in the system. When this superfluid state 
is disordered by proliferating both vortex species, the 
Coulomb phase obtains, and the domain wall excitations 
become particles identified with the gapped monopoles. 



III. VORTEX PHENOMENOLOGY AND 
PHYSICAL PICTURES 

The two-chargon gauge theory is amenable to a precise 
duality treatment, which gives a dynamical theory of vor- 
tices. A formal derivation is given in Appendix ^ Here, 
we proceed gently and consider a semiclassical descrip- 
tion of vortices in the superfluid. This discussion provides 
a physical picture that underlies the more formal treat- 
ments of subsequent sections. We argue that, depending 
on the microscopies, it is possible to have either a unique 
physical vortex or two distinct vortex species. (A more 
careful discussion requires the full dynamical theory.) As 
we vary the parameters in the superfluid phase, there is 
a transition occurring inside the core of vortex lines from 
a unique quantum core state to two degenerate but dis- 
tinct states. The "transition line" cuts the superfluid 
phase into two regions as indicated in Fig. |31 

It should be emphasized that there is no transition in 
the ground state of the system. The transition occurs in 
the properties of an inflnite vortex line, which is a very 
high energy excitation. But this has a bearing on which 
insulating state is obtained when vortices condense: If we 
have a unique vortex which then condenses, we obtain the 
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conventional insulator. On the other hand, as we argue 
below and in subsequent sections, if we have two distinct 
vortex species and both condense, we obtain the fraction- 
alized Coulomb phase (see FigO. It is also possible to 
have two vortex species but condense only one of them, 
and this will give an insulating state with broken discrete 
symmetry (bi <-> 62 symmetry in the two-chargon model 
or more physically the sublattice interchange symmetry 
in the microscopic boson model); in the simple lattice 
model Eq. {T)) this phase is probably not realized. 

We begin with an effective description of the superfluid 
phase, which is obtained by expanding the cosines in the 
classical action Eq. (TJ; the terms read schematically 

(V01 - af + {V(l>2 - af + (V X . 

Shifting a a' ~ a V(/>2, the field a' is massive, 
and at low energies we are left with one Goldstone mode 
[V(0i — <p2)]'^, as expected in the superfluid phase. 

To study vortices, we apply an external electromag- 
netic gauge potential that couples to the physical boson 
number density; this amounts to the replacements 

{V(t)i -a- ^A^^t? + (V02 -a+ ^A^^tf + (V x a)^ . 

Factors ±1/2 correspond to the ±qb/2 charges carried by 
the chargons. An hc/qi, vortex is obtained by requiring 

on going around the vortex. At large distances from the 
core we have 

V01 a ^Ae^t = 0, V02 - a + ^A^^t = . (12) 

Consider configurations with 0i^2 winding by 27rmi^2 re- 
spectively, f V0i^2 • dl = 27rmi_2- From the above equa- 
tions we conclude 

(j) Aext ■ dl = 27r(mi - 1112) (13) 
<j> a - dl = 7r(mi -(- TO2) • (14) 

Thus, there are different ways to realize a physical vor- 
tex in terms of the chargon fields. For a unit vortex, the 
relevant realizations are (mi,m2) = (IjO) and (0,-1) 
since these have the lowest core energy. Notice that 
the two realizations carry the gauge field flux = +7r 
and — TT respectively. More generally, different realiza- 
tions of a given vortex are characterized by values of 
the gauge flux $a that differ by multiples of 2ti. In- 
dependent of specific realizations, we note an important 
distinction between odd-strength vortices, which carry 
flux $a = Trmod (27r), and even-strength vortices, which 
carry = Omod (27r). 

In the above, the gauge field a is treated as nonconi- 
pact. On a formal level, the multiplicity of the physical 



vortex appears to be an artifact of ignoring compactness 
of the gauge field. Indeed, by moving a monopole across 
the system, which is a dynamical degree of freedom, we 
can change the gauge field flux by 2tt. For the preced- 
ing discussion of a single vortex, configurations (1, 0) and 
(0, 1) that differ in by 27r can mix with each other and 
produce a unique physical vortex (the mixing occurs on 
the level of local line segments). In this case, the only 
gauge flux distinction that remains is between even and 
odd vortices. 

However, it is crucial to note that the possibility of 
the mixing discussed above is a dynamical question. We 
argue below that it effectively does not occur near the 
Coulomb phase boundary. 

Consider a single physical vortex. As shown above, the 
vortex carries internal gauge flux of -I-tt or — vr depending 
on whether it is realized as a vortex in (pi or (j)2. Now 
consider a "domain wall" inside the core of such a vor- 
tex where the internal gauge flux changes from -I-tt to 
— TT. As the internal gauge flux changes by 27r at such 
a domain wall, we may identify it with a monopole con- 
figuration of the gauge field. In the Coulomb phase this 
monopole costs a finite non-zero energy. The monopole 
gap is expected to stay on approaching the transition to 
the superfluid. Consequently, close to this phase bound- 
ary but in the superfluid side, there will be a finite non- 
zero energy cost associated with the domain walls inside 
the vortex core. This immediately implies the existence 
of two distinct vortex species with the same physical vor- 
ticity which are distinguished by the sign of the internal 
gauge flux in the core. 

Upon moving away from the superfluid-Coulomb 
phase boundary, the domain wall gap in the vortex core 
will eventually close. The domain walls will then prolifer- 
ate inside the core leading to the destruction of the vortex 
core order so that there will then be a unique vortex for 
any given physical vorticity. 

What is the physical meaning of the two vortex 
species? To explore this, first note that the two vortex 
species may be described as being either a full 2tt vortex 
in bi (but not in 62) or vice versa. Thus the amplitude 
of bi is suppressed in the core of the 0i vortex while that 
of 62 is not, and the reverse happens for the (j)2 vortex. 
Thus, the order in the core of the vortex corresponds 
to spontaneously breaking the symmetry of interchange 
of 61 and 62- But what is this symmetry in terms of a 
more microscopic description in terms of the underlying 
bosons? The answer clearly depends on the model con- 
sidered. For the microscopic boson Hubbard model of 
Sec.m the two chargon species 61.2 are associated with 
sites that belong to the two distinct sublattices of the cu- 
bic lattice on whose links the physical bosons live. Thus, 
the chargon interchange symmetry corresponds simply to 
interchange of the A and B sublattices in the underlying 
model. We therefore conclude that this sublattice inter- 
change symmetry (present in the microscopic model) is 
spontaneously broken in the core of a vortex line in the 
superfluid phase close to the transition to the Coulomb 
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insulator. The two different vortex species simply cor- 
respond to the two possible degenerate ordered patterns 
for the core state. 

An almost similar situation also obtains for bosons at 
half- filling (see Section 0) on a cubic lattice in a super- 
fluid state close to the transition to a Coulomb insula- 
tor. There the two different vortex species again have 
broken sublattice symmetry in the core which may be 
simply identified with checkerboard density wave order. 
This kind of order in the vortex core is very natural for 
bosons at half-filling — thus a Coulomb insulator is per- 
haps a more natural possibility at half-filling once the 
superfluidity is destroyed. 

Having understood the vortex structure in the super- 
fluid phase, we can now turn around and address the 
question of how to view the Coulomb phase in terms of 
the excitations of the superfluid. To destroy the superflu- 
idity we need to proliferate vortex lines. With two vortex 
species present, first consider proliferating just one kind 
of vortex. The resulting insulator will then inherit the or- 
der in the core of the proliferated vortex species. In other 
words, the resulting insulator will break a symmetry (for 
instance, the sublattice interchange in the boson Hub- 
bard model of Sec^J. On the other hand, we may pro- 
liferate both vortex species with equal amplitudes. The 
resulting insulator will then not inherit the broken sym- 
metry of cither vortex core. This is the Coulomb phase. 
The monopoles are clearly the remnants of the domain 
walls between the two ordering patterns along the vortex 
core (see Fig. The photon may then be associated 
with fluctuations of a 'roton' formed by combining to- 
gether a loop of one vortex species with a loop of an 
antivortex of the other vortex species. To understand 
this, we note that in the gauge theoretic description of 
the superfluid phase such a roton has no net vorticity 
but carries internal gauge flux of 27r. In other words, it 
represents a tube of magnetic flux. Once the two vortex 
species proliferate, the fluctuations of the roton may be 
described in terms of fluctuating magnetic field lines, and 
this physically describes the photon. 

The above physical pictures are elaborated in precise 
formal terms in the following sections. 



IV. DUAL VORTEX DESCRIPTION 

To make the above discussion of the vortex core struc- 
ture and its implications for proximate insulating phases 
more precise, it will clearly be useful to have a dual for- 
mulation which focuses on the vortex line degrees of free- 
dom. For the issues at the focus of the present paper — 
namely a dual description of the Coulomb phase — it is 
necessary to have a dual formulation that clearly brings 
out the presence of the two distinct vortex species with 
ordered cores. This is most conveniently obtained by du- 
alizing the gauge theoretic description in terms of the two 
chargon fields. Such a formulation may be obtained by an 
extension of the duality transformation for the usual lat- 



tice XY model in (3-|-l)D. This model describes bosons 
at some large integer filling on a three dimensional lat- 
tice. Since even this duality is possibly unfamiliar to 
most readers, we first review it below and show how var- 
ious familiar Mott insulating states of bosons at integer 
filling may be described in terms of vortex line conden- 
sates. We will then generalize this duality to incorporate 
the non-trivial core structure of the vortex lines when 
there are two vortex species, and show how the Coulomb 
phase results from their proliferation. 



A. Review of vortex description for (3+l)D XY 
model 

Here we review duality for generic 3D quantum bosons 
at integer filling. The dualities are performed for the 
Euclidean action, which in this case is (3+l)D XY model: 



SxY[<^)^] = -/3^cos(V^<; 



(15) 



To avoid clutter, we take equal spatial and temporal 
couplings. The formalism is standard^SiiS but we want 
to br ing ou t its physical interpretation as was done in 
Refs. [nil 3 for (2-t-l)D. Our treatment below makes ex- 
plicit the correspondence of the dual variables to vortici- 
ties, and enables us to obtain a dual perspective on many 
familiar Mott insulating states. The primary purpose of 
this review is to provide a background for the subsequent 
duality treatments of more complicated cases. 

We start with the Villain form for the partition func- 
tion 



bi J —TT 



(16) 



We work on a (3+l)D hypercubic lattice; space-time sites 
are labeled by lower-case Latin letters (?,...), and the 
lattice directions are labeled by Greek letters (/i, ...). 
Integer-valued fields pifj_ ensure 2tt periodicity in the an- 
gles. For each configuration in the statistical sum, we 
define vorticities 



(17) 



In (3-|-l)D, g^i^ describe vortex worldsheets. Here and 
below, the lattice derivatives are indicated only schemat- 
ically, but the precise meaning is readily recovered in each 
case. 

As usual for 0(2) models, we want to separate the 
spin wave and vortex parts. This is achieved by dividing 
configurations {pifj,} into classes with the same vorticity. 
Two configurations {pi^} and {p^^} belong to the same 
class if they can be related via p-^ = pi^ + V^iV with an 
integer- valued field TV^. Using the latter, we now extend 
the (p integrals over the whole real line and obtain 
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Zv = 



^ nOO nOQ 

/ [Djzf.] (^iV^jp) exp 

1 J — OO 



(18) 
(19) 



r 



Here, we introduced a real- valued field j^^, which can be 
interpreted as particle current; is any member of 

the class with given vorticity g^^^. The prime on the 
sum over the vorticities qip,i, indicates that these satisfy 
schematically d{qfj^udXf^ A dx^) — 0, which reads as an 
integer-valued constraint 



1 _ 



(20) 



for each direct lattice cube or equivalently each dual lat- 
tice link. Here and in what follows, we specifically work 
in (3+l)D and make use of the fully antisymmetric tensor 
EpCT/ii^ to connect between direct and dual lattice objects; 
summation over repeated indicies is implied unless spec- 
ified otherwise. The above constraint means that there 
are no sources for vortices. We can also specify the vor- 
tex worldsheets by an integer- valued field Fjpa- living on 
the dual plackets (upper-case letters specify dual 

lattice sites): 



1 



"^aFpa = 



(21) 

(22) 



The current conservation constraint ^ pjfi = can be 
solved by introducing a rank-2 antisymmetric field Qipa 
living on the dual plackets 



(23) 



2 27r 

with factor 27r introduced for convenience. We now have 



(24) 



I,p<a 



SO vorticity appears explicitly in the statistical sum. At 
this stage, we can integrate out the field Qp^ and obtain 
a description in terms of the vortex worldsheets with 
specific long-ranged interactions, and consider different 
phases of the vortex system as determined by this in- 
teraction. We expect that modifying the interactions at 
short distances will not change the physics of the vari- 
ous phases but only positions of the phase boundaries, 
and we therefore consider generalized models with added 
local "fugacity" terms 



'S'fug. — — 



1 

2A 



1 

2A 



E 

I,p<(T 



I pa 



(25) 



In order to describe the vortex system in more famil- 
iar terms, we instead retain the gpo- field and proceed 
as follows. The integer-valued constraint of no vorticity 
sources written in terms of Fjp„ , Eq. (|22|l , is handled by 
introducing a J7(l) variable c/p on each dual lattice link 
in the manner 



S{V^Fp^) cx 



/ dCpCXp [iCpVc-Fpa] 
J — 7r 



(26) 



(no sum over p). We can now perform the summation 
over vorticities obtaining 



Zv[X] = / [D0ip,]' / [i?c/p]exp 



2^E-^'v 



X exp 



A ^ cos(V 

I,P<(T 



(7^p 



(27) 



where the cosine stands for the appropriate Villain form 
and is given in Eq. H23|l. The only approximation in 
the above duality transformation (besides Villain-izing) 
is the added vortex fugacity. The resulting dual theory 
has a compact U{1) gauge field coupled to a noncom- 
pact rank-2 field. This is a generalization of the (2-|-l)D 
duality, where the dual theory has a U{\) scalar field 
coupled to a noncompact gauge field. The vorticity g^^ 
of the original angle variables is precisely the integer- 
valued electromagnetic field tensor Fp^ of the compact 
electrodynamics in the c/p variables. In a Hamiltonian 
formulation, vortex lines are identified with the electric 
field lines, e.g., = F^x = qyz- Thus, e*'^"'" is conju- 
gate to Eux and can be viewed as a vortex line segment 
creation operator. 

For later convenience, we also exhibit the dual lattice 
action using 'soft-spin' vortex fields: 

^vort = E [*Jp*?+P,^*/+*:P*/- e-'S- +C.C.] 

I,p<<J 

+ ^V^(|*/pn + «5](ep.p.V,0p,)^ (28) 



ip 



Here 5'}^ 



jicjp creates a vortex line segment on the 
dual lattice link; y(|*P) = s|*|2 -|- ^41^-14 + ug|^|6 + 
• • • controls amplitude fiuctuations of the field From 
the preceding derivation, the rank-2 field g/po- describes 
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boson density fluctuations via Eq. H23() . The first term in 
the action represents vortex Hue "hopping" by extending 
across a placket. The boson density acts as a source 
for the rank-2 gauge potential seen by the vortices (this 
gauge potential produces the familiar Berry phase when 
a vortex line moves in the superfluid). 

The behavior of the compact U{1) gauge theory with- 
out the rank-2 field is well understood, and we can use 
this as a starting point to develop intuition about the vor- 
tex system. This line of thinking can be interpreted as 
considering screened vortices first. Indeed, for a charged 
boson system we would conclude that the field 2p(t is mas- 
sive and can be ignored. To obtain a faithful description 
of neutral bosons, e.g. to correctly reproduce the low- 
energy modes, we need to include the "gauging" by the 
rank-2 field. This can be accomplished by a semiclassical 
analysis just as in the Ginzburg-Landau theory. 

We now summarize the dual description of the familiar 
phases of 3D quantum bosons at integer filling. 

Superfluid phase of the bosons is identified with con- 
finement in e"^'' . Vortex excitations in the superfluid are 
the gapped electric fleld lines. The rank-2 fleld gpo- is 
essentially free, and its single propagating mode (which 
can be found e.g. by doing quantum mechanics for this 
field) is precisely the phonon mode of the superfluid. 

Mott Insulator phase is obtained by proliferating 
vortices, so e'^'^p is "condensed", or more appropriately 
"deconfined". The rank-2 field Qpa is gapped out, and in 
the process "eats" the photon modes (an analog of the 
Anderson-Higgs phenomenon), so there are no gapless ex- 
citations, as expected. Gapped charged excitations of the 
Mott insulator are represented by monopoles of the dual 
gauge field, and the discreteness of charge is encoded in 
the corresponding quantization condition for the rank-2 
field, which is a generalization of the flux quantization in 
the Ginzburg-Landau theory. 

Indeed, let us consider a monopole worldline. Fix a 
time slice (of our Euclidean path integral) and evaluate 
the boson number 

ntot = J jrd^r = ^ j Cpr{V^Qyz + V^fl^a: + "^zQxy) ■ 

(29) 

In the far field away from the monopole we have e.g. 

Qxy = ^xCy ~ VyCx = Bz , (30) 

so the right hand side coincides with the outgoing mag- 
netic flux divided by 27r, which is equal to the en- 
closed monopole number. Thus, monopole worldlines in 
the dual description are precisely worldlines of integer- 
quantized charges in the direct description. 

Z2 Practionalized Insulator: We can readily gen- 
eralize the above description of the Mott insulator to 
discuss more exotic insulator phases. Z2 fractionalized 
phase is obtained as a condensate of doubled vortices 
while single vortices are not condensed. The dual action 



generically allows double- vortex hopping terms like 

-A2 ^ C0s[2{y pCa -V aCp~ Ql pa)] ■ (31) 

/,p<(T 

When A is small while A2 is large, we expect e*^'^f to be 
deconfined but not e*'^''. Topological excitations in this 
phase are monopoles of the field e*^'^'' , and these are seen 
to represent charged excitations carrying fraction 1/2 of 
the boson charge. Odd-strength vortex line excitations 
correspond to vison line excitations of the Z2 phase. 

B. Dual formulation with two vortex species 

We now show how to incorporate the physics of two 
distinct vortex species with ordered cores in the dual 
formulation. As mentioned above, this is most conve- 
niently done by dualizing the gauge theoretic descrip- 
tion in terms of two chargon fields. The duality trans- 
formation is a straightforward but tedious extension of 
the methods of the previous subsection. We therefore 
relegate the technical details to Appendix ^ and focus 
here on the result and the physics contained in it. As 
expected on the basis of the semiclassical discussion in 
Sec, mil the analysis in Appendix ^yields a dual theory 
with two vortex fields '^^^^ and ^'^^^ corresponding to 
the two possible ordering patterns in the core. Further, 
these vortices are allowed to turn from one flavor into the 
other, and the corresponding domain wall appears as an 
additional field T in the dual theory. 

A generic 'soft-spin' action for the two vortex fields can 
be expressed in the form 

fc=l,2 

- A„,5](T;T,+p*(;)*vI;(2)+c.c.) 

+ Y.U{\TA^). (32) 

In the two-chargon theory, both ^!\^ and ^ create 
a physical vortex line segment, but it is realized as a 
vortex in bi in the first case and an antivortex in 62 in 
the second case, see Eq. (I13|l . The vortex interaction 
terms V and W need not be specified at this stage but 
are important when determining the vortex condensate. 
Am represents quantum tunneling between the two core 
states and is written here as the domain wall hopping, 
while U represents the domain wall energy cost. 

Two broad possibilities need to be considered in the 
superfluid. The flrst possibility is that domain walls are 
energetically cheap and proliferate inside vortex cores. In 
this case the two vortex fields lock to each other. In the 
integer filling case considered here, this produces a unique 
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physical vortex, and the proper low-energy description is 
then in terms of a single vortex field (the situation at half- 
integer filling is more complicated, see the next section). 

The second possibility is that domain walls are ener- 
getically costly and remain gapped inside vortex cores. In 
this case the field T has only short ranged correlations 
and may be integrated out. The two vortex fields are 
proper low-energy degrees of freedom. We can now ask 
what happens when these vortices condense, e.g. when 
^' = becomes an unstable state of the potential V. De- 
pending on the interaction W, the condensate either has 
just one of the two vortex species, or has both vortex 
species proliferated simultaneously. The former case oc- 
curs when W > 0, and we write schematically (^f^^^) ^ 0, 
(xI/(2)) = 0. The resulting Mott insulator then inherits 
the order in the core of the vortex. There are no 
gapless modes: as discussed in the previous subsection, 
the rank-2 field Higgses out the gapless modes of the 
single-species worldsheet condensate. In the microscopic 
boson Hubbard model on the lattice of corner-sharing 
octahedra, this Mott insulator breaks A ^ B sublattice 
interchange symmetry. 

On the other hand, for the interaction W of the oppo- 
site sign, it is favorable for both vortex species to prolifer- 
ate simultaneously, schematically = (^f^^^) ^ 0. In 
this case we obtain fractionalized Coulomb insulator with 
a gapless photon, a gapped monopole, and two species of 
gapped chargons. The gapless photon obtains since the 
rank-2 field Qpa can Higgs out only the \I'(i)*v[/(2)* 'part' 
of the two-species worldsheet condensate. The vl/(i)*vl/(2) 
part does not couple to gpo-; it represents string objects 
with short-range interactions only — indeed, it is pre- 
cisely the roton formed by combining a vortex of one 
kind with an antivortex of the other kind. To show that 
its fluctuations lead to a gapless photon mode, we focus 
on just the fluctuations of the 'phases' of the 
two vortex species: 



(33) 



[Note that compared with Appendix 1X1 we use 

c(i), '^#'(2) = -c(2)^ while the monopole field 6'(") is the 

'phase' of the domain wall particle ~ e^^' .] Expand- 
ing to quadratic order in small gradients of 'rf^^\^'-^\ 
the resulting 'elastic' action has the following schematic 
structure: 



S = (V X 



fl)2 + {V x'^ 



(34) 
(35) 



Here g represents the rank-2 gauge field. The ellipses 
refer to the kinetic energy term for the g and various an- 
harmonic corrections to the above quadratic action. It 
is now easy to see that the combination ^^^^ -I- "^^^^ is 
rendered massive while the combination ^^^^ — 'rf^^^ re- 
mains massless. The latter precisely represents a linear 
dispersing gapless photon with two transverse polariza- 
tions in three spatial dimensions. Thus as promised the 



photon is indeed associated with fluctuations of the roton 
formed from a vortex of one kind and antivortex of the 
other kind. The gapped T field simply corresponds to 
the monopole excitation of the Coulomb phase. Finally, 
the "monopoles" of the proliferated dual fields 5"^^^ and 
^ff^' are the two chargons. This completes our formal 
description. We also point that Appendix IdI develops a 
similar description of the Ising P* phase in two spatial 
dimensions that incorporates the physics of two distinct 
Ising domain walls. 

We conclude by emphasizing that the vortex prolifer- 
ation transition is likely first order, in which case the 
lattice action for vortices does not have a formal contin- 
uum limit, but this does not affect the presented physical 
picture of the phases. 



V. BOSONS AT HALF-FILLING 

We now consider bosons at half-filling in some detail. 
We have in mind some generic Hamiltonian 



Ub 



E 



(36) 



The physical picture presented in the Introduction and 
elaborated in the previous section carries over readily and 
is perhaps more appealing in this case. In particular, the 
ordering pattern in the vortex core when there are two 
species simply corresponds to checkerboard density order 
of the bosons. This is a very natural "insulating" core 
for the vortices to develop when the bosons are at half- 
filling. Formally, this core structure may readily be seen 
to arise in the same semiclassical treatment as in Sec. lIIII 
We first break up the boson field into two 'chargon' fields 
as before. This leads to a description in terms of a com- 
pact U{1) gauge theory coupled to two chargon fields. A 
physical vortex in the superfiuid order then corresponds 
to having a full vortex in bi (but none in 62) or a full 
antivortex in 62 (and no vorticity in hi). Thus, there 
are still two types of vortices which convert into each 
other at the locations of monopoles of the gauge the- 
ory. Following the closely related discussion in Ref. [3 in 
two dimensions, it is readily seen that the two vortices 
correspond to the two possible patterns of checkerboard 
density wave ordering in the core. (The precise reason 
is simply that the two vortices differ in the sign of the 
internal gauge flux. For bosons at half-filling, this gauge 
flux couples linearly to the difference of the boson densi- 
ties on the two sublattices; hence there is checkerboard 
density wave ordering in the core). 

As for bosons at integer filling, the insight into the vor- 
tex structure can be given a precise form by employing 
a duality transformation. The dual action is derived in 
AppendixEl The Coulomb phase is again obtained when 
both vortex species condense with equal amplitude. The 
domain walls in the ordered vortex core survive as gapped 
monopoles, and the photon emerges exactly as before (as 
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the fluctuations of proliferated rotons made by combin- 
ing a vortex and an antivortex of two different species). 
The gapped chargons are again "monopoles" of the two 
vortex fields that appear in the dual theory. Thus the 
main structure of the description of the Coulomb phase 
is unaffected from that appropriate at integer filling. 

The primary differences arise when we consider the 
nature of the possible confined phases. These may be 
obtained from the Coulomb phase by condensing the 
gapped monopoles. Such a condensation leads to confine- 
ment of the gapped chargons and the disappearance of 
the photon from the spectrum. However, as can be antic- 
ipated from experience in two dimensions, the monopoles 
now transform non-trivially under lattice symmetries. 
This leads to broken lattice symmetry in the confined 
phases. In particular, various bond or box density wave 
phases are possible. These will be explored directly in 
the next section. 

Here we specify more precisely the modifications to the 
two-vortex action Eq. (|32() due to half-filling on a cubic 
lattice and sketch some implications for the vortex core 
structure. As discussed in Appendix lO there is an ad- 
ditional static vector potential X^^ that couples to the 
domain wall hopping term — the term in Eq. (|32|l — 
see Eq. (|C4p . This vector potential encodes the monopole 
Berry phases due to half-integer boson density. It is spec- 
ified by the corresponding fluxes = VpX^ — V^-^p 
(modulo 27r) through the faces of the dual cubes. The 
fluxes are nonzero on the spatial plackets only, and equal 
to 7r/3 when oriented from the A to B sublattice of the 
direct cubic lattice as shown in Fig. El The derivation of 
this result is given in the next section, and further details 
can be found in Appendix O 

As discussed above, the two flelds '^^^^ and corre- 
spond to the two checkerboard charge orders inside vor- 
tex core. We can now ask about the effect on the core 
dynamics due to the static vector potential X'^ seen by 
the domain walls. The speciflc details of the domain wall 
motion such as change in the dispersion are not imme- 
diately relevant for the low-energy physics if the domain 
wall remains gapped. Indeed this is why the modifica- 
tions induced by half-filling do not affect much the de- 
scription of the Coulomb phase itself. But they become 
more important when the monopole gap closes at a phase 
transition from the Coulomb to a confined phase, or when 
the domain wall gap closes inside vortex core. 

Deep in the superfluid phase when the vortex line is 
stiff, we expect a single minimum in the domain wall dis- 
persion, and when the walls proliferate inside the core we 
obtain a unique vortex. However, for a vortex line that 
fluctuates sufficiently and explores the three-dimensional 
space so that the effects of X^ are felt, the bottom 
of the domain wall band can split. This is the situa- 
tion in the Coulomb phase, where domain walls become 
monopole particles, and one finds two low-energy propa- 
gating monopole modes (see the next section for details) . 

Consider now what happens when we have two domain 
wall band minima inside the vortex core in the super- 



fiuid phase, and the domain wall gap goes to zero. One 
possibility motivated by our analysis in the next section 
is that domain wall proliferation leads to a VBS order 
inside the core. In this case, while the fields 5'^^-' and 
vl'(^) are locked, there can be several degenerate ways for 
this to happen corresponding to the number of favorable 
VBS states inside the core (we need to be even more 
careful here since e.g. for the columnar order the favor- 
able states likely depend on the vortex line orientation). 
This would mean that our choice of the basic fields 
and ^E"^^^ is not a very good one, and the analysis needs 
to be reconsidered, perhaps introducing different vortex 
fields. When a single such vortex condenses, the corre- 
sponding VBS state results, but it may be also possible 
to condense several such vortices simultaneously possibly 
leading to new exotic states. Even more exotic possibility 
is for the domain walls to become critical inside the core, 
which might occur due to the one-dimensionality of the 
setting. In the present work, we have only touched upon 
these possibilities and have not pursued any systematic 
studies. 

We now comment on the possibility of a direct anal- 
ysis of a more 'microscopic' vortex theory for bosons at 
half-filling and the connection with the above two- vortex 
description. Indeed such a direct analysis is very useful 
in two spatial dimensions The bare vortex action at 
half-filling is described in Appendix IbI As in (2-|-l)D, 
the half-filling manifests itself by "frustrating" the vor- 
tex propagation with additional static Berry phases. In 
(3-|-l)D, this is encoded by placing half of a magnetic 
monopole inside each dual cube in the compact gauge 
field part that describes the vortex motion. As men- 
tioned in Appendix at present we do not know how 
to analytically treat such frustrated gauge theories and 
how to connect directly with the two-vortex description 
like it was done in Ref. 14 for (2-f 1)D. Our earlier discus- 
sion and the analysis in Appendix [CI provide an alterna- 
tive route in two dimensions'^ that extends reasonably 
to three dimensions, and we expect the resulting action 
with two vortex species to describe the vortex physics at 
low energies in some regime of parameters. (A limited 
but direct attack is possible on the bare vortex action as 
well^) 



VI. VALENCE BOND SOLIDS OF BOSONS AT 
HALF-FILLING ON A 3D CUBIC LATTICE 

In this section, we show how valence bond solid 
phases emerge due to monopole condensation out of the 
Coulomb phase at half-filling. The crucial ingredient is 
the monopole dynamics produced by the static vector 
potential X'^. To that end it is convenient to specialize 
to the hard-core limit in which case we may view the 
bosons as representing spin-1/2 moments living on the 
cubic lattice. In this limit, the static vector potential 
may be viewed as arising from Berry phases present in 
a path integral description of a spin-1/2 moment. Be- 
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cause of the importance of this result, which is an exten- 
sion to three dimensions of Haldane's calculation in two 
dimensions we give a direct derivation in the context 
of spin-1/2 systems. It is in this context that the va- 
lence bond solids were analyzed in the work of Read and 
SachdevfSi providing us with familiar grounds. 

A. Monopole action and Berry phases 

Our starting point is the 'modified' Sachdev- 
JalaberliiSiS^'S^ lattice model for spin-1/2 system on a 
cubic lattice in the vicinity of the collinear Neel state. 
This is defined in terms of a bosonic field Zia which car- 
ries spin-1/2 and lives on the sites (denoted «, j, ) of a 

(3-|-l)D space-time lattice (a =t,i is a spin index). We 
will refer to the quanta of this field as spinous. They are 
the analog (for the spin system) of the chargon fields in- 
troduced earlier for boson models. The spinous are min- 
imally coupled to a compact U{1) gauge field a^^. The 
action for the Sachdev-Jalabert model reads 




X exp {-S^ - Sa - Sb) , 



S. = -/3^[z^e"-z,+^-Kc.c.]-f ••• , 
Sb = i^^ViO'iT ■ 

i 

The compactness of the gauge field is encapsulated by 
summing over integer- valued fields u^^y. The last term 
represents a Berry phase that encapsulates the spin-1/2 
nature of the moment at each site. The rji appearing 
in this term is +1 on the A sublattice and —1 on the 
B sublattice of the spatial cubic lattice. Connection to 
microscopic spin models and mot ivation may be found 
in the original references l23l24l2R With appropriate de- 
formations of Sz, this model can also describe bosonic 
systems at half-fnteger filhng (see Refs.EHl). 

Here we are interested in paramagnetic states with 
gapped spinous (interpreted as chargons in bosonic sys- 
tems). In such cases we may integrate out the spinon 
fields completely. This effectively amounts to dropping 
the Sz term from the action. (As we are interested in the 
general structure of the paramagnetic phases, we will not 
worry much about effects such as generating other short 
ranged interactions for the gauge fields). The remaining 
action is that of a compact U{1) gauge theory in (3+l)D 
on a cubic lattice but in the presence of the Berry phase 
term. The latter has a simple physical interpretation in 
a Hamiltonian language. It simply corresponds to the 
presence of a static background gauge charge of strength 



±1 on the sites of the spatial cubic lattice. The sign 
alternates from one sublattice to the other. Indeed pre- 
cisely the same kind of gauge theory arises in studies of 
two-dimensional spin-1/2 quantum paramagnets. 

In the (2-t-l)D case, it is well-known that the alter- 
nating static background charge leads to broken transla- 
tional symmetry in confined paramagnetic phases. Our 
goal is to extend that analysis to three dimensions. 
To that end we perform a duality transformation on 
Sa + Sb- In the absence of the background charge, this 
is just the familiar electric-magnetic duality. The com- 
pact t/(l) gauge theory then becomes a theory of point- 
like monopoles that are coupled to the dual fluctuating 
non-compact U{1) gauge field. (The non-compactness of 
the dual gauge field is because we have thrown out the 
spinous — which are the gauge 'electric' charges — in the 
original theory. Indeed retaining them is equivalent to re- 
taining monopoles in the dual gauge field). The presence 
of the background charge means that there is dual mag- 
netic fiux emanating out of the center of each cube on 
the dual lattice on which the monopoles live. This flux 
alternates in sign from one cube to the next and frus- 
trates the monopole hopping. We therefore expect the 
duality transformation on Sa + Sb to yield a theory of 
monopoles with frustrated hopping that is coupled to the 
dual non-compact gauge field. This is explicitly shown 
in what follows. 

We proceed by writing 

/oo 
-OO 

X exp[i ^ fi^i^CVf^ajy - V^^a^ - 27rui^,y)] . 

ffj,i, can be interpreted as the electromagnetic field tensor. 
We classify configurations {wi^i/} by their monopole four- 
currents 

Two configurations {u^^} and {m^j^} belong to the same 
class if they can be related u'^^ = Ufj,^+\7 ^V^ — \7 ^V^ with 
an integer- valued field V^. This allows separation of the 
Gaussian and topological defect parts of the statistical 
sum, since the integration now extends over the full 
real line. 

To treat the Berry phase term, wc find a static /^^ 
satisfying 

V^Spr = V,/°, , (38) 

f'. ^ ~^iSf.rV.m~5,rV^m) . (39) 

(As will be clear below, the proportionality constant 
uqK/ P is in fact fixed by the above conditions to be 
1/12; the specific writing is chosen for convenience when 
reusing the present analysis in Appendix lUl) The first 
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condition is possible since = r]i5f,,r has zero diver- 
gence. It allows us to write the Berry phase term as 



(40) 



lJ.<v 



We can now bring out the topological defect part by 
changing variables to f^i, — — /^^. In terms of the 
variables /^^ the action has the same form as in the ab- 
sence of the original Berry phases except for additional 
contributions 



a)h° 




Si 



-E 

K ^ 



/m^C + *E/"-2™^- ■ (41) 



Due to condition (|38l) the 'phase' term Sb depends on the 
monopole configuration only. The second condition H39() 



allows us to write the term S'„q as 



•^"0 - E 



no 



(42) 



Here = Viy f^v can be viewed as the gauge charge 
current. This follows by integrating over the field; 
Appendices 1X1 and O exhibit this explicitly for U{1) mat- 
ter fields, while in the absence of any matter we get 
^uffiu = 0. Thus, Sno can be viewed as a static staggered 
chemical potential seen by the z fields, and as long as it 
is much smaller than the spinon gap, it can be ignored 
for the low-energy properties. 

We now exhibit that solves the two conditions 
Eqs. lISHIl and llSni). Only the electric field parts E° = f°^, 
Ey, and E'^ arc nonzero and satisfy 



(43) 
(44) 



This fixes E'^ uniquely to have value 1/6 on each bond 
from the A to B sublattice (the value no/f3 is also fixed 
at this stage). 

The phase term Sb can be expressed in terms of the 
monopole currents Jp™'' by first transforming from f^^ 
to i)~ 



per- 



f,0 

J 2 f^'^P'^ 27r ' 



(45) 



The latter has spatial components f)° y , l}^^ , [}°^ only, 
which can be interpreted as fluxes through the faces of 
the dual lattice cubes. The total outwards flux is +2n 
for a dual cube surrounding an A-sublattice site of the 
direct lattice and — 27r surrounding a i3-sublattice site, 
and in each case the flux is divided equally among all six 
faces. This is shown in Fig. Since the total flux out 
of each cube is ±27r, we can write 



(1° 



VpX^ 



V„X° mod 2tt 



(46) 




even 



FIG. 6: Description of the monopole Berry phases for bosons 
at half-filling on a simple cubic lattice in terms of a static 
gauge potential seen by the monopoles hopping on the dual 
lattice, a) fj^a gives a flux of ±7r/3 through each spatial 
placket of the dual lattice with the flux oriented from the 
A to B sublattice of the direct lattice, b) Our gauge choice 



for e'^p that realizes the fluxes e^^'p'^ (Eq.0BJ. 



Xp is a static field on the spatial links of the dual lattice. 



and one choice is shown in Fig. 03. We finally obtain 



p 



(47) 



and interpret this as a monopole Berry phase. This is ap- 
propriate since as we have explained the remaining con- 
tributions to the low-energy action have the same form 
as in the absence of the original Berry phase Sb- 

We can offer the following geometric representation of 
the monopole Berry phase. Consider a monopole space- 
time worldline and trace out the spatial path traveled by 
the monopole. The Berry phase is given by the flux of f)" 
(see Fig. through any surface pulled onto this closed 
path. This is (3-|-l)D generalization of Haldane's (2+l)D 
result^ for the monopole Berry phase. In (2-f 1)D, the 
Haldane's result leads to destructive interference among 
monopole worldpaths unless monopoles are quadrupled. 
In (3-|-l)D, the interference pattern encoded in f)^^. is 
more subtle, and requires us to study the monopole hop- 
ping in the static vector potential X^; this is detailed 
below. 

We now complete the duality mapping for the model 
Sa + Sb without any matter field in order to exhibit the 
relevant structure of the monopole action; the relation- 
ship of this to the original Sachdev-Jalabert model will 
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be commented upon later. 

Using the standard result for the monopole action in 
compact QED we obtain 

Here, Lp is a noncompact dual gauge field such that 

//it/ = i-^Lvpa^ pi^a I i^T^)- The prime on the sum in- 
dicates that the monopole currents satisfy continuity 



VpJp™'' = 0. Integrating over Lp we would obtain an 
action in terms of the monopole worldlines only. Instead 
of doing this, we transform the monopole action into a 
more familiar form by considering a generalized model 
with added monopole fugacity term 



1 



'fug. 



2A„ 



(m)N2 
P > 



(48) 



The constraint V pJp [I) = on each dual site / can 
be solved by introducing a U{1) field 9^/^\ and upon 
summing over Jp™"* the final result reads 



/OO /"TT 
[DLjp]' / [D6j 
-OO J —IT 



(™)lg-S[ip,e('"'] 



S 



A„ ^ cos(Vp0(") - X° - Lp) . (50) 



e'" can be interpreted as a monopole creation opera- 
tor; the field 6'^™-' is coupled to the noncompact (dual) 
gauge field Lp, and there is also static "frustration" 
coming from the Berry phases for monopoles in our spin- 
1/2 system. 

The relationship of the above action to the original 
Sachdev-Jalabert model is as follows. The action was 
derived in the absence of the spinon field z (e.g. taking 
the spinon gap to infinity) . Now, z carries electric charge 
of the original gauge field a^, and as usual in the electro- 
magnetic duality, it acts as a magnetic charge for the dual 
gauge field Lp. Thus, to have a complete correspondence 
with the original spin model, we need to allow monopoles 
in the dual gauge field Lp, and these monopoles need to 
be spin-1/2 particles. Ignoring the spinon field makes the 
dual gauge field noncompact as in the above action. 



B. Monopole condensation patterns 

We now study possible phases of the above monopole 
action assuming throughout that the spinous are gapped. 

Clearly, for small Am the monopole field is gapped. In 
this case, the dual gauge field is free, and we obtain a 



L'^(l) spin liquid phase of the original spin-1/2 problem. 
This phase has gapped spinous, gapped monopoles, and 
a gapless photon, and respects all lattice symmetries. 

To analyze possible monopole condensates, we closely 
follow Ref. lUl. We first ignore fluctuations of the dual 
gauge field Lp and study a frustrated XY model defined 
by the static vector potential What we are essen- 
tially doing here is averaging over fast fluctuations im- 
posed by Once the structure of the resulting slow 
fluctuation description becomes clear, we will restore the 
held Lp. 

Continuous-time soft-spin action for the frustrated XY 
model is 



R 



(RR' 



R 



V{Wr?) 
(51) 

with some potential y(|$fip) = rol^j? P + uo|<i'i?|^ H ■ 

We work on the discrete spatial lattice labeled by R (dual 
to the original spin system lattice), since the lattice is 
crucial at this stage. 

The frustration is encoded in the monopole hopping 
amplitudes = te^-^RR' , corresponding to the fluxes 
\f' on the plackets as shown in Fig.EJi. Our gauge choice 
is shown in Fig. and details as follows: 



(49) tR,R+x 




tR,R+V 
tR,R+z 



1 



Diagonalizing the kinetic energy, we flnd two low-energy 
modes with normalized real-space wave functions 



*i(i?) 

*2(i?) 



/2(3 - V6) 
1 - ( \/3 - \/2)e'' 
2(3 - V6) 



V2 



(52) 
(53) 



Thus, we find two monopole excitations carrying different 
lattice momenta. 

At the kinetic energy level, any linear combination 



$(i?) = (/.i«'i(i?) + 02«'2(i?) 



(54) 



is at the bottom of the monopole band, and there is a 
continuum of states for the monopoles to condense to. 
Non-linear terms will lift the degeneracy. This can be an- 
alyzed near monopole condensation transition by treating 
01 and (f)2 as slowly varying flelds and deriving Ginzburg- 
Landau theory of these. 

By examining the action of the lattice symmetries, the 
resulting Ginzburg-Landau functional is required to be 
invariant under the following transformations (in our spe- 
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cific gauge): 
Ty 

R90°,Rxy 
R90°, Rxz 



for EiT' = + C.C., etc.; while the 

monopole current is 



^1 1 



V2 ' 

; 



(^2 1 



V2 



V2 



The 90° rotations are about the lattice points on which 
monopoles reside. The hopping Hamiltonian is invariant 
(modulo gauge transformations) under the above lattice 
symmetries if we also complex-conjugate tRR' i*RR' be- 
cause of the "staggered flux" pattern — this is the reason 
for having 0*s on the right hand sides. There is no di- 
rect significance to the seeming "time-reversal breaking" 
in the specific gauge; the physical states of the spin sys- 
tem that we study here (with gapped spinous) are time- 
reversal invariant. 

The simplest invariants are \(f>i\'^ + |'/'2p and 

4 (i</>in^2p + \ci>i\^m') - 6i</>ih02i' - {r.hr 

The latter can be given a more clear form 

/8(0i, </'2) = N^N^ + N^Nl + NlNl , 

with 



(55) 



(56) 



We write down continuum action for the two- 
component complex field </>(i?, t) that respects the above 
symmetries 



5= 



low 



drSR 



\{V p-iLp)cj>\''+U{\(i>\'')+vsh{ct,) 



-St 



(57) 

with some potential [/(|</>p) = r\(f)\^ -\- U4\(f)\'^ -\- ■ ■ ■ . We 
have also restored the dual gauge field Lp and included 
some generic kinetic energy 5*^ for Lp. 

We focus on confining paramagnets that obtain start- 
ing from the Coulomb phase and condensing single 
monopoles so that {<p) ^ 0; this happens when r < 
in the potential. [It is also possible to have paired con- 
densates such that {(f>^<7(f)) ^ but (0) = in some 
parameter regime — such states correspond to U{1) spin 
liquids with broken translational symmetry; these are not 
considered in any detail here.] Ground states are selected 
by minimizing vsls{(f>), and the sign of wg determines the 
character of the resulting phase. 

Each state can be characterized by the expectation val- 
ues Nx, Ny, defined above. For example, the spatial 
monopole density in a given state is given by 

|$(i?)P = + -1= [{-irN.. + {-l)yNy + i-lYN^] {f>8) 



the monopole kinetic energy is 



E. 



(m) 



(59) 



J-. 

ff") = l{t 



(60) 



for Ji'"-' = i{tR^R+-t^*j^^R+x - c.c). 

We also quote transformation properties of TV = 
{Nx: Ny, Nz) under the lattice symmetries: 



Tx : 


{Nx,Ny,N,) - 




Ty: 


iNx,Ny,N,) - 


-> {Nx,-Ny,N,) 


: 


{Nx,Ny,N,)- 


-> {Nx,Ny,--N,) 


R^Q" ,Rxy ■ 


{Nx,Ny,N,)- 


{Ny,Nx,N,) ; 


R9Q°,Rxz ■ 


iNx,Ny,N,) - 


^ iN,,Ny,Nx) . 



We remind that the rotations are about points R of the 
lattice on which monopoles reside, which is dual to the 
original spin model lattice. Given the above information, 
we can characterize broken lattice symmetries when N 
obtains an expectation value. 



1. «8 > 

In this case, there are six ground states 

(7V,,Ar^,iV,) = (±1,0,0), (0,±1,0), (0,0,±1). 

(61) 

In a given state, the monopole density is the same in 
every other plane perpendicular to a fixed lattice axis. 
For example, the state {N^, Ny, N^) = (0, 0, 1) has an in- 
creased density on the 2;-even planes and decreased den- 
sity on the z-odd planes, and is illustrated in Fig. UJi. 
In the original spin model, these dual planes are crossed 
by direct lattice bonds, and on a crude level there is an 
increased bond energy crossing the z-even planes. Thus, 
this state corresponds to a columnar valence bond solid 
with dimers oriented in the z direction and in the par- 
ticular registry. This identification is also supported by 
an analysis of the dual gauge field fluxes induced by 
the monopole currents and the connection between these 
fluxes and the bond order parameternS^ The above six 
states correspond to six possible columnar states on the 
cubic lattice. We expect that this is the phase realized 
when monopoles condense in the specific lattice model 
Eq. (|50|l . based on the analysis of the lattice XY model 
and the large Am limit of the action Eq. H50|) . 



2. vs<0 



In this case, there are eight ground states 
(iV„iV„7V,) = -^(±l,±l,±l) 



(62) 



(each of the three signs can be chosen independently). 
These states have monopole density oscillating in all 
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FIG. 7: a) Schematic picture of columnar VBS state obtained 
when Us > 0. We have an increased monopole density (and 
also energy density) on the shaded z-even planes, which is in- 
terpreted as having valence bonds of the original spin model 
preferentially crossing these planes. Arrows along the placket 
edges indicate monopole currents, b) 3D box VBS state is 
obtained when us < 0. For clarity, only the planes with in- 
creased monopole density are shown (a:;, j/, z all even), and the 
dimers resonate around the cubes centered where these planes 
meet. 



three directions and correspond to three-dimensional 
"box" valence bond sohds analogous to box states in 
two dimensions. For example, the state {N^, Ny, N^) = 
-^(1,1,1) has maximal monopole density for x, y and 
z all even, and is illustrated in Fig. [Tja. In the original 
spin model, this state has dimers resonating around di- 
rect lattice cubes surrounding these dual lattice points. 
The above eight states correspond to eight possible ways 
to register such 3D box states. 



Discussion 



The above considerations lead us to the following pic- 
ture. In the spin- 1/2 system we expect VBS order pa- 
rameter to be a three- vector 



^(-1)^5, • Sr- 



^VBS 



drd^ 



We can write down a Ginzburg-Landau functional for 4" 

(V^*)2-fr*2+u(#2)2+^,(^4_^^4^^4) 

(64) 

where besides 0(3)-invariant terms we have also included 
quartic interaction with cubic anisotropy w, since our sys- 
tem resides on the cubic lattice. Thus, the VBS phases 
in three dimensions are described by an effective 0(3) 
model with cubic anisotropy. Depending on the sign of 
the anisotropy, the three-vector 4" can either point along 
one of the axes [e.g. ^E" (x (0, 0, 1)], which leads to colum- 
nar order with 6 degenerate states, or along a diagonal 
[e.g. ^ cx (1, 1, l)/\/3], which leads to box order with 8 
degenerate states. 

The connection with the preceding analysis is provided 
by the identification 

N = 4)^34, . 



VBS 



(65) 



(63) 



Indeed, 5* has the same transformation properties un- 
der lattice symmetries as N . Thus, the three-vector ^ 
is written with the help of the two-component complex 
field 4> precisely as in the CP^ representation of the 0(3) 
nonlinear sigma model. As we know, a proper descrip- 
tion of the 0(3) model in the CP^ language contains be- 
sides the spinor field </> also a compact C/(l) gauge field, 
and the monopoles of this gauge field are conventionally 
identified with the hedgehog configurations of the 0(3) 
vector. This CP^ gauge field is precisely our gauge field 
L in the 0- field action Eq. (|57|l . and as discussed earlier, 
for a faithful description of the original spin- 1/2 system 
we need to allow spin-1/2 carrying magnetic monopoles 
in the gauge field L, which correspond to the gapped 
spinous. 

The identification is complete when we realize, extend- 
ing the discussion in a recent preprint by M. Levin and 
T. Senthilj^e that the VBS order parameter field ^vbs 
indeed has hedgehogs which are precisely the spinous! 
This is dictated by the microscopic lattice origin of 
^VBSi which needs to be brought back into consideration 
(this microscopic detail is absent in the continuum action 
SvBS, Ea. l64|l . Consider for example the columnar phase 
and construct a hedgehog in the field TV by putting to- 
gether the corresponding domains of columnar dimer or- 
der. More precisely, we are constructing a "hedgehog" 
of the corresponding six-state discrete model. When we 
join these domains, there is an unpaired site left which we 
identify with a spinon residing on one sublattice; and if 
we construct an antihedgehog, we find a spinon residing 
on the other sublattice. Of course, we pay huge domain 
wall energies consistent with the fact that the spinous are 
confined in the VBS phase. 



D. Connection with the vortex description 

We conclude by describing how to obtain the 
above picture directly from the two-vortex action 
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Eqs. and (|C4|I at half-filling. We use hard-spin fields 
and notation as in Appendices lAICI Consider phases 
with gapped chargons obtained by proliferating both vor- 
tex fields c^^^ and c'-^-' . Ignoring the chargons completely, 
we can expand the cosines containing V x c, and after 
simple analysis arrive at the monopole action Eq. H5()|l 
with Lp = c'p + cjf' ■ The real utility of the two- vortex 
action is that it allows us to incorporate the chargons on 
equal footing in the dual description as "monopoles" in 
the vortex fields c*-^-* and c*-^-*. In particular, it formalizes 
for the easy-plane case the purported picture that we are 
to include monopoles in the gauge field Lp to obtain the 
full physical description; no such formalism is known for 
the Heisenberg case. 

VII. CONCLUSIONS 

In this paper, we described a physical mechanism for 
the appearance of the fractionalized Coulomb phase in 
bosonic models. This was accomplished by addressing 
the question of how to view this phase as a vortex con- 
densate emanating out of a proximate superfluid phase. 
We were led to suggest one route whereby there are two 
degenerate but distinct physical vortices in the superfluid 
phase, and the superfluid is destroyed by proliferating 
both species of vortices with equal amplitude. The two 
vortex species correspond to vortices whose cores (which 
are one-dimensional systems) are ordered and break a 
discrete symmetry. A simple example is provided by 
bosons at half-filling where the vortex core develops insu- 
lating checkerboard charge density order. The two vor- 
tex species correspond to the two checkerboard states 
in the core. Domain walls in the core order survive as 
gapped excitations when the two vortex species prolifer- 
ate equally. We argued that these domain walls may be 
identified with the monopole excitations of the Coulomb 
phase. Roton loops formed by combining a vortex loop 
with one core order with an antivortex loop with the 
other core order were shown to correspond to loops of 
the emergent magnetic flux. The fluctuations of these 
loops correspond to the gapless photon excitation of the 
Coulomb phase. 

In general, the microscopies that produces degenerate 
but distinct physical vortices is model-specific. But quite 
broadly, we can view these vortices as having ordered 
cores, where the order is that of some nearby insulat- 
ing phase with a broken discrete symmetry. For generic 
bosons at half-filling, one simple candidate is the checker- 
board charge order, while for the specific corner-sharing 
octahedra model discussed in Sec. |n] the candidate has 
checkerboard 'octahedron' order. 

The primary analytical tool used in the paper to put 
meat into this description is a duality mapping to vor- 
tex variables which is an extension of the familiar boson- 
vortex duality to (3+l)D. We extensively developed such 



a description for the quantum systems of interest and 
used it to discuss the physics. To discuss the Coulomb 
phases, we showed how the presence of two vortex species 
can be readily incorporated into the dual theory. The ex- 
citation spectrum of the Coulomb phase was then derived 
in this dual formulation which provides detailed confir- 
mation of the physical pictures. This theory also allows 
us to identify and describe all basic nearby insulating 
states in a unified setting. 

For bosons at half- filling (closely related to spin- 1/2 
quantum antiferromagncts), apart from the fractional- 
ized phases, various confining phases are possible which 
break lattice translation symmetry. The latter were 
analyzed by extension to (3-|-l)D of the methods of 
Haldancr^ Read, and Sachdev^^ (originally developed for 
lower dimension). This allowed us to discuss possible va- 
lence bond orders in this system, and their relation to 
the fractionalized Coulomb insulating state. 

In the present study of bosonic superfluids and prox- 
imate insulators at half-filling, we concentrated on vor- 
tices with charge-ordered cores. It is conceivable instead 
to have valence bond order inside the core in some regime 
of parameters in the superfluid phase, and it would be in- 
teresting to develop a broader picture of the core dynam- 
ics and its consequences. The above considerations pro- 
vide further explicit examples of nontrivial effects arising 
from the vortex core physics for disordering transitions of 
XY ordered systems, echoing similar phenomena in two- 
dimensional systems." -^'^^ In particular, we showed that 
there can be transformations occurring inside the core, 
and that proliferating ordered cores can produce unusual 
phases. The importance of the core physics may have 
broader implications for experimental studies of strongly 
correlated systems. 
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APPENDIX A: DUALITY TRANSFORMATION 
FOR TWO-CHARGON GAUGE THEORY EQ. 

Here we perform duality analysis of the two-chargon 
gauge theory Eq. (TJl. Our final result is the two- vortex 
action Eq. (|32|l at integer filling. We first identify topo- 
logical defects in the model. Villain-ized partition func- 
tion for the model is written as 
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exp{i^ [ji ■ (V01 a 27rpi) + j2 • (V02 - a - 27rp2)] + i ^ fitiui^^iCLy - V^a^ - 2™^^^,)} . (Al) 



r 



Real fields ji , j2 , /^^i/ can be interpreted as the matter 
currents and the electromagnetic field tensor. 

We divide configurations {pi,p2,u^y} into classes 
equivalent under integer-valued transformations 



p'i = Pb + VNt,-V, 6=1,2, 



(A2) 



with integer fields Ni, N2, V. These classes are specified 
by two vorticities 



and monopole currents 



t(t) 



The latter satisfy continuity 



Vp = 



(A3) 



(A4) 



(A5) 



while the vorticities satisfy dq*-''-' — du, which can be 
written as an integer-valued constraint 



1 



(A6) 



The meaning of the last equation is that monopoles act 
as sources and sinks for both vorticities simultaneously. 

For each allowed configuration of vorticities and 
monopole currents, we can now perform the summation 
over all {pi,P2,M^y} in the corresponding class using 
Eqs. (|A2p . This effectively extends the integration vari- 
ables (j)i,4)2, dp, over the whole real line, and the integrals 
give the delta function conditions 

[^(V • ji)][5(V • j2)][<5(V,/p, - j'lp - j2p)] (X 



In the last line, we changed variables to = ji ± j2. 
Note that not all conditions in the first line are indepen- 
dent, and the precise meaning of the delta functions is 
given by the second line. Thus, we can completely elim- 
inate j + . 

The topological defects enter through 



«7r^(j+(pi +P2) +j (Pl -P2))+i'2TT^fpuUf,^ 

p<v 



p<i' 



We now solve the constraint V • j =0 with a rank-2 
field g7p£r defined on the dual plackets 



^pupa 



(A7) 



It is also convenient to pass from the variables /^^ to 
their dual g+g. via 



Ipu = 7:e 



pupa 



(A8) 



J"*" is obtained from g"*" by an expression similar to 
Eq. (|A7|) . Finally, it is convenient to specify vortici- 
ties g^t' by the corresponding integer- valued fields , 
6=1,2, 



^ pa n'^papv'Jpu 



V7 p{b) _ T(m) 
^ c-f^pcr ~ "^p 1 



(A9) 



just as we did for the (3-Hl)D XY model in Eq. if^ . 
Putting everything together, the partition sum reads 



Z 



E 



L 1 p<7 • 1 per ' 1 p J 



I pa 



I pa I 



e 



+„(F<i) + f(^))p„+B-^(F 



^AlO) 



We can now in principle integrate out the fields and with sources and sinks on the monopole worldlines. The 
Qpg., and obtain a theory in terms of vortex worldsheets 
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vortex interactions are such that combinations f 
enter with short-ranged interactions similar to screened 
vortices, while F'^^^ — F^^^^ enter with long-ranged inter- 
actions of unscreened vortices. This is expected since the 
chargon field combination + (f)2 is gauged while (jji — 02 
is gauge neutral. We expect the physics to remain un- 
changed upon adding local vortex and monopole fugacity 
terms 



Just as in the case with the XY model in Sec. IIV Al 
to get a better intuition about the theory at hand, we 



instead introduce two compact U{1) gauge fields 



„(2) 



„(1) 



and a f7(l) scalar field dj' that implement the con- 
straints Eqs. (|A9IA5p and are the appropriate conjugate 
variables. 



'S'fug. 



1 

2A 



I,p<cr 



1 



2\r. 



Ip 



(All) 



The final expression reads 



/TT poo 
-TT J — OO 

X exp[A„, J2 cos(Vp0(") - c« - cfj) + A ^ cos(Vpc(i) - V^c^ - 0+ , - g," J 
Ip l,p<'y 

+A ^ cos(Vp42) - V,c(2) _ 0+ ^ + J] . (Af 2) 

I,p<a 

I 



As usual, cosines stand for the appropriate Villain forms. 
At this stage, the field is massive and can be inte- 
grated out. The variables of the resulting dual theory are: 

i) two compact U{1) gauge fields e*"^'^' and e*'^*^', which 
can be viewed as vortex line segment creation operators; 

ii) monopole matter field e*^*"' that carries both dual 
gauge charges; iii) noncompact rank-2 field q~ that fur- 
ther "gauges" the vortex fields c*^^^ and c'^'^'^ and describes 
the physical boson density fluctuations. Equation l|32() in 
Sec. IIV fJl displavs a 'soft-spin' version of this dual theory 
by writing 

vl/W*^e'4^.', vl/(2)*^e-'^^'.', TJ^e^^^"' (A13) 

(we also dropped the superscript on the rank-2 gauge 
field g~ that survives in the low energy theory). 

The analysis of the possible phases from this dual per- 
spective can be done as in the XY case in Sec. IIV Al 
by first considering a theory without the q~ field (i.e., 
"screened vortex" theory which would obtain if the two 
chargons were also coupled to a fluctuating external elec- 
tromagnetic field.) The field q~ is then included semi- 
classically and provides correct count of the low-energy 
modes. 

Thus, the screened vortex theory has a phase for 
A 1, A,„ ^ I in which both fields c'-^-' and c'-^-' are 
deconfincd while the monopole field 6''^™^ is gapped. The 
screened vortex theory has two photons in this phase. 
Including fluctuations in the rank-2 field g^, it obtains 
a mass and eats in the process one photon, so there re- 
mains precisely one photon [V x (c^^-* -I- c'-^-')]^. In the 
physical boson model, this is the fractionalized Coulomb 
phase at the focus of this paper. 



On the other hand, in the regime A ^ 1, A„i ^ 1 
monopoles are also proliferated and kill the photon, so 
this becomes the conventional Mott insulator of bosons. 

Finally, for small A ^ 1 both dual gauge fields c^^^ and 
c*^2) are confined. The dual electric field lines which are 
the vortices of the physical boson model are gapped, and 
this corresponds to the superfluid phase of bosons. As 
we vary the parameter A,„ inside this phase, there is a 
transition in the properties of these line excitations from 
a unique core state to two degenerate core states. 

The dual formulation reveals some details of this core 
transition. Specifically, we consider a single vortex in the 
screened vortex theory — this allows us to avoid compli- 
cations of coupling the vortex dynamics to phonons and 
focus on the core physics. 

It is convenient to use Hamiltonian formulation 

ifscrcenod = U,^ Y.'^R" E COs(V0(™) - C^^) - C^^)) 

+u^^{e\ + el) - i« E [cos(V x c^^^) + cos(V x c^^^) , 

with the Hilbert space constraints 

V ■ei = V ■e2 = nR . (A14) 

Here u,„ and are bare monopole gap and hopping 
amplitude, while is bare vortex core energy per unit 
length and is vortex hopping amplitude. The super- 
fluid phase of the boson model corresponds to the the 
confining phase u^j ^ ty of i?scrconcd- The electric field 
e = ei — 62 is the conserved physical vorticity, and we 
study a strength- 1 line of e oriented in say x direction. 
In the limit ^ ty, we have a straight line e^, = 1, but 
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its segments can be realized as either {eix,e2x) = (1,0) 
or (0,-1). Each juncture (0,-1) (1,0) has a +1 
monopole and (1,0) (0,-1) has a —1 monopole, and 
the quantum dynamics is determined by the competition 
between the u„i and t„i terms. Writing erf — +1 for (1, 0) 
segment and af = —1 for (0, —1) segment, we obtain the 
following ID Hamiltonian for the straight line 



-Hcore ^ Z^'^l ^l+l ^-^7.0"; 

/ I 



(A15) 



This is simply the one-dimensional quantum Ising model, 
and has two phases: For Um < tm-, obtain a paramag- 
net in ct^ and the ground state is unique, schematically 
CT^ = 1. In this case, the monopoles proliferate along the 
line, and the e = 1 line is unique. On the other hand, for 
Um > tm, 'Vfe obtain a ferromagnet with two degenerate 
ground states = +1 and cr^ = —1. The monopoles 
remain gapped along the line, and we have two distinct 
e = 1 lines. The core Hamiltonian can be also written 
in dual Ising variables, and this gives a description in 
terms of domain wall particles which we identify with 
the monopoles. 

The above analysis is valid in the limit Uy S> i« , but we 
expect the two regimes and the transition between them 
to extend throughout the confined phase of i?scrocncd- In 
the superfluid phase of the original boson model, we con- 
jecture that the corresponding vortex core transition line 
extends all the way to the insulating phases as shown in 
Fig. El 



APPENDIX B: DUALITY ANALYSIS FOR 
BOSONS AT HALF-INTEGER FILLING 



Here we generalize the analysis in Sec. IIV Al to bosons 
at half-integer filling and study the resulting dual theory 
directly. At half-filling, we simply make the replacement 



,2 

JlT 



(Bl) 



in the final expression Eq. (|27|l while the rest remains un- 
changed. As in (2+l)D, noninteger boson density enters 
as an external "field" seen by the vortices, except that 
vortices are now lines and the gauge potential is rank-2 
tensor. 

A convenient formulation is obtained by finding static 



such that 



(B2) 



If flp^ are viewed as fluxes through the dual lattice plack- 
ets, 0° has total flux of tt coming out of each spatial cube. 



Shifting by g° we get 

/OO /'IT " 2^ 



X exp 



A ^ cos(VpCCT — V, 

I,p<a 



(T^p Ql pa 



^ pa 



Thus, there is a static contribution to the rank-2 gauge 
potential corresponding to placing half of a magnetic 
monopole inside each spatial cube. When a vortex line 
evolves in space-time, this gives a Berry phase contribu- 
tion of TT for each cube in the 3D volume swept by the 
line. 

As in our discussion for bosons at integer filling, it is 
helpful to first consider the compact U{1) part with only 
the static 0°^. This "frustrated" lattice gauge theory 
has not been studied to our knowledge. At present, we 
do not have direct tools for approaching this problem 
similar to the ones used for frustrated systems with global 
symmetries. We can still describe what we expect to 
happen in such a model. 

Two phases of the frustrated lattice gauge theory are 
clearly identified. For small A, we expect confinement in 
gzc/p^ which corresponds (upon restoring fluctuations in 
0pcr) to the superfluid phase of the bosons. 

For large A, we expect deconfinement in e'^'^''', but we 
need to be more specific since there can be different pat- 
terns of deconfinement in the presence of 0pg.. For the 
'hard-spin' action written above, we can proceed by con- 
sidering first the classical ground states and obtain the 
following picture. For each spatial cube, the total out- 
going flux of e*^^'''^''"^'''^''"^^''^ must be tt modulo 27r 
because of the compactness of the gauge field. For an 
individual cube, the lowest energy is obtained by divid- 
ing this flux equally among the six faces. On the lat- 
tice, we put such cubes with the outgoing placket fluxes 
of -|-7r/6 and cubes with the outgoing fluxes of — 7r/6 in 
a checkerboard pattern, and there are two degenerate 
ground states. This structure survives as a phase in the 
frustrated lattice gauge theory. For the original bosonic 
system, this phase corresponds to checkerboard charge 
density wave. Indeed, when we allow fluctuations in Qp^, 
the above staggered character is imprinted on the fluxes 
0pcr, thus producing staggered boson density. The two 
states correspond to two ways to register this CDW on 
the lattice. 

Other deconfinement patterns are also possible, but 
likely require some deformation of the above simple ac- 
tion to be stabilized. For example, it is possible to have 
a deconfined state with an energy density wave, and this 
would correspond to some VBS phase for the original 
bosons. 

Of direct interest here is the possibility of deconfine- 
ment in e^'^'p that produces two photons in the low-energy 
description and does not break any lattice symmetry. Al- 
lowing fiuctuations in gp^, only one photon remains, and 
this would give the sought for fractionalized Coulomb 
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phase at half-integer fiUing. We know that this can in- 
deed happen in the original bosonic system, and Appen- 
dices H and |0 provide an indirect route to describe such 
a phase. It would be interesting to explore this in the 
above frustrated lattice gauge theory in more detail and 
develop direct analytical tools for treating such theories. 



APPENDIX C: DUALITY FOR TWO-CHARGON 
THEORY AT HALF-INTEGER FILLING 

In this Appendix, we derive dual description for a 
two-chargon gauge theory corresponding to bosons at 
half-integer filling. Our analysis is a direct extension to 
(3-hl)D of Appendix E of Ref. E which should be con- 
sulted for further details. Following Ref. 0, we arrive 
at a path integral which has the same form as at integer 
filling Eq. (|A1|) but with replacements 

{jlirf {jur - mVif , {j2irf {32ir - nol]if , 

and an additional term in the action 

Sb =i'^rjiair ■ (CI) 

i 

In the limit K ^ (i.e. in the absence of the gauge field 
dynamics), any choice of the "chemical potential" hq for 
the on-site gauge charge will give the same result due to 
constraint jur + j2iT = const. This is no longer true for 
finite K, but the qualitative behavior is expected to be 
insensitive to no, and we will use this freedom to make 
the structure of the dual theory more transparent. 

We proceed as in Appendix IXI while treating the phase 
term Sb as in our analysis of the monopole Berry phase 
for spin-1/2 systems in Sec. I VII Specifically, we take f^^ 
defined by conditions Eq. H38I39|) and make the corre- 
sponding change of variables /^i^ = /^^ -I- . In terms 

of the variables ji^, j2/j, f^i^, we obtain the same action 
as at integer filling except for additional contributions 

In particular, upon integrating out 01, 02, a^i, we have 
j+ — Vuf^iij- The first two terms cancel each other if 
Kuq/P =1/12 (see Eg. I42|l . The importance of the spe- 
cific form for /^^ is that for general uq the first two terms 

add up to ^ J2i Vijtr^ which can be viewed as a stag- 
gered chemical potential seen by the gauge charges. As 
long as this potential is small, it effectively averages out 
and does not modify the low-energy modes; in the fol- 
lowing, we set (5no = 0. 

The remaining phase term can be expressed entirely in 
terms of the monopole currents Jp™' 




FIG. 8: Quantum Ising model that realizes P* phase in two 
dimensions. Spins are located on the link-centered sites of a 
square lattice. Each shaded diamond indicates the ?7-term in 
the Hamiltonian Eq. ljDl|l . 

with defined in Section IVII This term needs to be 
added to the dual action expressed in terms of the topo- 
logical defects Eq. (|A10|) . Thus, we have managed to 
absorb all effects of half-integer filling into the monopole 
Berry phases. To obtain the final action in terms of the 
dual fields 0*^™), c^^\ c*^^', we simply take the expression 
at integer filling Eq. (|A12p and make the replacement 

Vpe*'™) ^ Vp6^"'^ - a:)] . (C4) 

APPENDIX D: ISING P* PHASE IN (2+l)D 

Toy model: Here we consider a toy model which is 
a 2D quantum Ising version of the 3D boson model of 
Sec.m As we will see, some physical aspects are similar 
in the two problems, while the formalism is much simpler 
in the Ising case and gives some guidance for the bosons. 
Ising spins reside on links of a square lattice, and the 
Hamiltonian is 

H = -jY^alat,-TY^af-UY^X{at. (Dl) 

{W) I r ler 

The lattice is shown in Fig. |S1 The spins can be viewed 
as forming a network of corner-sharing diamonds, and we 
introduced four-spin interaction U on each diamond. For 
simplicity, we take U > 0. 

The phase diagram of the model can be analyzed as 
in Sec. ^1 When the Ising coupling dominates, J ^ 
r, U, the system is ferromagnetically ordered. When the 
transverse fields dominate, T,U ^ J, the system is a 
conventional Ising paramagnet. 

The intermediate regime U ^ J ^ \/TU is most inter- 
esting and realizes an unusual topologically ordered para- 
magnet (Ising P* phase). The excitation spectrum in this 



23 



phase contains a gapped Ising vortex (vison) and gapped 
Ising matter fields. The argument is similar for the Ising 
and boson models. In the large U limit, the ground state 
sector is determined by the projection Yiier '^i ~ ^ '^'^'^ 
all r, and the effective Hamiltonian in this sector is 

I □ 

with K = J'^ /U . This is just the familiar Ising gauge the- 
ory. For K '^T the theory is confining and corresponds 
to the conventional paramagnet, but for X > F we have 
a deconfined phase and this is the Ising P* phase. 

The topology of the phase diagram is the same as in 
Fig. ^ with the correspondences superfluid — > ferromag- 
net, conventional insulator conventional paramagnet, 
and fractionalized insulator P* paramagnet. We are 
led to ask similar questions on the relationships among 
the phases, in particular, how to view the P* phase com- 
ing from the ferromagnet. In the ferromagnet, excitations 
are domain walls, and we can think of a paramagnet as 
a result of proliferating these. The P* phase features the 
vison, which is a topological point excitation, and we ask 
how this appears in the domain wall condensate. 

To push the analogy somewhat further, in (3-|-l)D the 
appearance of the Coulomb phase with the monopole 
excitation is 'natural' when an 0(3) spin model is dis- 
ordered without proliferating hedgehogs (gapped hedge- 
hogs then become the monopoles). As discussed in the 
main text, the appearance of the Coulomb phase in a 
nominally 0{2) model is more subtle in terms of the vor- 
tex line defects of the ordered phase. Similar situation 
occurs in (2+l)D, where the appearance of the Z2 frac- 
tionalized phase with the vison excitation is natural when 
an 0(2) model is disordered by proliferating double vor- 
tices but not single vortices, but the appearance of such 
Z2 phase is not so clear in a nominally Ising model. We 
argue that the key to this puzzle lies in the presence of 
two degenerate but distinct domain walls in the ferro- 
magnetic phase near the P* phase. 

Direct analysis of domain walls: In the present 
model, we can perform a direct analysis in the limit of 
small r, i.e., in the upper right corner of the phase dia- 
gram Fig. ^ 

When r = in the microscopic model, we find that for 
each square lattice placket the corresponding Z2 'flux' 
<I>P = o'J'2C23f''34C''4i is Conserved. This simplifies the 
analysis, since in each sector with fixed $'s the Hamil- 
tonian decouples into two quantum Ising models. In- 
deed, let us work in the basis and fix a sector {$p}. 
Let us find a specific realization {cr^^?^} of the fluxes, 
<i>P = o'J'j^^ 0-23° Vll" VlJ''^ Wrr'} bclougs to this scctor 
if we can write cr^^, = S^S^,a^l^\ and we can label all 
states in the sector by the new Ising variables {S^}. The 



action of the Hamiltonian in this sector reads 
= - E Jaa'S^S^, - f/E -51 

{AA') A 

- E jbb'S%s%, -uy^^Sg . 

{BB') B 

Here A and B refer to the two sublattices of the square 
lattice. Ising couplings Jaa' and Jbb' depend on the 
sector {$p}; for example, Jaa> is determined by con- 
sidering the square lattice placket [1^,2, 3^', 4] that has 

A A' as its diagonal: Jaa' = •^(''■12°^ '''23°^ + ^''14°'' '^43°^)- 

For arbitrary J/U (with F = 0) the ground state sector 
has no fluxes, <i>p = 1, and Jaa' — Jbb' = 2 J. As 
we change J/U, we have a simultaneous Ising ordering 
transition for spins Sa and Sb- The gap to the nearest 
sector is at least of order J in the ferromagnetic phase 
and of order J'^ /U in the spin-disordered phase. 

Consider the ferromagnetic phase. In the sector with 
no fluxes, excitations are domain walls for the Ising vari- 
ables Sa or Sb- When translated to the original variables 
domain wall in Sa is a physical domain wall that 
passes through the B sublattice sites, while a domain 
wall in Sb is a physical domain wall that passes through 
the A sites. In a sector with two visons, one finds that 
the lowest energy state has an A-type domain wall (de- 
fined by frustrated bonds J aa' S\S\,) connecting the two 
visons, and also a -B-type domain wall connecting the vi- 
sons. Thus, visons act as sources for both domain walls 
and are linearly confined. 

Consider now the spin-disordered phase. It can be 
viewed as a condensate of both domain walls, and has 
two Ising matter fields Sa and Sb with the mass of or- 
der U . The different sectors are accounted for as vison 
excitations with the vison gap of order J'^ /U . 

For nonzero F, the different sectors will mix, but the 
above description of the phases will remain true as long 
as F is much smaller than the gap to the nearest sector. 
Small nonzero F will induce energy-energy coupling be- 
tween the two Ising models; this will affect the nature of 
the transition, either making it first-order or driving to a 
different universality class; however, the above physical 
picture with two distinct domain walls remains. 

Phenomenological gauge theory description sim- 
ilar to the two-chargon theory is also possible. Again, it 
relies on the observation that in the P* phase we find two 
distinct Ising matter excitations. The phenomenological 
Hamiltonian is 

~r E '^rr' - -^E ^12*^23^34 <^41 ; 

{rr') □ 

it contains two Ising matter fields Tir,T2r coupled to an 
Ising gauge field arr' ■ The Hilbert space of the theory is 
defined by 

Tl^^2. n "rr' = 1 • (D3) 
r'Gr 
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We emphasize that this model is not derived from the 
Hamiltonian Eq. (|D1|I (in particular, the reader should 
not be confused by our reuse of the letters for the coupling 
constants). Rather, this model is to be considered as 
capturing the relevant physics of the P* phase. 

The analysis of the phenomcnological gauge theory 
proceeds exactly as in the two-chargon theory, and the 
phase diagram has the same topology as in Fig. [31 The 
ordered phase has ferromagnetic order in the 'physical' 
spin = rf^T|^. In the ferromagnetic phase, we have do- 
main wall excitations, and a domain wall segment cross- 
ing a link t^t^i = — 1 can be realized as either a domain 
wall in Ti {a^^iTl^T^j., = —1) or a domain wall in T2. Z2 
flux (vison) acts as a source for both domain walls, and 
there is a linearly confining potential between a pair of 
visons in the ferromagnetic phase. On a physical domain 
wall, however, vison costs only finite energy and can hop 
along the domain wall; this hopping represents quantum 
tunneling between the two microscopic realizations of the 
hopped segment. Vison can be thought of as a point 'do- 
main wall' for the order along the line which is the Ising 
domain wall of the ferromagnetic phase. 

We expect two possibilities: One is that visons prolif- 
erate along the domain wall, in which case the domain 
wall is unique. The other possibility is that visons remain 
gapped inside the domain wall, and there are two distinct 
domain walls. A precise formulation is to consider an ex- 
ternally imposed domain wall in the ferromagnetic phase. 
As we change the parameters in the ferromagnetic phase, 
there is a transition from a disordered domain wall with a 
unique ground state to an ordered domain wall with two 



degenerate ground states, and this domain wall ordering 
transition is of (l-l-l)D Ising type. When we have ordered 
domain walls and proliferate both with equal amplitudes, 
the P* phase results. 

The mathematical formalism for the above picture is 
transformation to dual variables, which is readily done in 
the Hamiltonian language. The dual Hamiltonian is 

i?2is,dual = ~^ '^^(^'1 2 ^2 3 ^34 ^41 + -^12-^23 -^34-^41) 

□ 

{RW) R 

the Hilbert space is defined by 

R'<£R R'eR 

The two domain wall variables are denoted as C, D, while 
the Z2 flux variable is denoted as V: Cf^jii = (J^j.iTi^tI^, , 

Drr' = <r'T-2rT|r': and = (rf^f^j ^f si^ ii- Compare 
this with the dual formulation Eq. (|A12|I of the two- 
chargon theory. The ferromagnetic phase in the dual 
language corresponds to confinement in both dual gauge 
fields C and D. In the P* phase, both these fields become 
deconfined while the dual matter field V where the two 
connect remains gapped. 
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